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, Abstract. We report on ten type I X-ray bursts originating from GX 17+2 in data obtained with the RXTE/PCA 

■ in 1996-2000. Three bursts were short in duration (~10s), whereas the others lasted for ~6-25min. All bursts 

showed spectral softening during their decay. There is no evidence for high-frequency (>100 Hz) oscillations at any 
phase of the bursts. We see no correlations of the burst properties with respect to the persistent X-ray spectral 
properties, suggesting that in GX 17+2 the properties of the bursts do not correlate with inferred mass accretion 
rate. The presence of short bursts in GX 17+2 (and similar bright X-ray sources) is not accounted for in the 
Q |- current X-ray bursts theories at the high mass accretion rates encountered in this source. 

We obtain satisfactory results if we model the burst emission with a black body, after subtraction of the persistent 
pre-burst emission. The two-component spectral model does not fit the total burst emission whenever there is 
a black-body component present in the persistent emission. We conclude that in those cases the black-body 
contribution from the persistent emission is also present during the burst. This implies that, contrary to previous 
suggestions, the burst emission does not arise from the same site as the persistent black-body emission. The 
black-body component of the persistent emission is consistent with being produced in an expanded boundary 
k>( i layer, as indicated by recent theoretical work. 

Five of the long bursts showed evidence of radius expansion of the neutron star photosphere (independent of the 
spectral analysis method used), presumably due to the burst luminosity reaching the Eddington value. When the 
burst luminosity is close to the Eddington value, slight deviations from pure black-body radiation are seen at 
energies below ~10keV. Similar deviations have been seen during (long) X-ray bursts from other sources; they 
can not be explained by spectral hardening models. 

The total persistent flux just before and after the radius expansion bursts is inferred to be up to a factor of 2 
higher than the net peak flux of the burst. If both the burst and persistent emission are radiated isotropically, 
this would imply that the persistent emission is up to a factor of 2 higher than the Eddington luminosity. This 
is unlikely and we suggest that the persistent luminosity is close to the Eddington luminosity and that the burst 
emission is (highly) anisotropic (£~2). Assuming that the net burst peak fluxes equal the Eddington limit, applying 
standard burst parameters (1.4 Mq neutron star, cosmic composition, electron scattering opacity appropriate for 
high temperatures), and taking into account gravitational redshift and spectral hardening, we derive a distance 
to GX 17+2 of ~8kpc, with an uncertainty of up to ~30%. 

Key words, accretion, accretion disks — binaries: close — stars: individual (GX 17+2) — stars: neutron — X-rays: 
bursts 
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X-ray bursts were discovered in 1975 from the source 
4U 1820-30 (Grindlay & Heise 1975; Grindlay et al. 1976). 
It was realized soon thereafter that these were thermo- 
nuclear runaway events on the surface of neutron stars 
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(Woosley & Taam 1976; Maraschi & Cavaliere 1977). 
Another kind of X-ray bursts was found (together with 
the above type of bursts) from MXB 1730—355 (later re- 
ferred to as the Rapid Burster), which were suggested to 
be due to accretion instabilities. The former and latter 
kind of bursts were then dubbed type I and type II, re- 
spectively (Hoffman et al. 1978a). 

The main characteristics of type I bursts (for a review 
see Lewin et al. 1993) are: sudden and short (~1 s) increase 
in the X-ray flux, exponential decay light curve, duration 
of the order of seconds to minutes, softening during the 
decay (attributed to cooling of the neutron star surface), 
(net) burst spectra reasonably well described by black- 
body emission from a compact object with ~10km radius 
and temperature of ~l-2keV, and total energies ranging 
from ~10 39 to 10 40 erg. When the luminosity during the 
burst reaches the Eddington limit (i.e., when the pressure 
force due to radiation balances the gravitational force), the 
neutron star photosphere expands. Since Lb oc R 2 T ff 4 , 
when the radius of the photosphere, R, expands, the ef- 
fective temperature, T e ff , drops, with the burst luminosity, 
Lb, being constant (modulo gravitational redshift effects 
with changing R) at the Eddington limit, LEdd- Bursts 
during their radius expansion/contraction phase are there- 
fore recognizable by an increase in the inferred radius 
with a simultaneous decrease in the observed tempera- 
ture, while the observed flux stays relatively constant. 

The emission from a (hot) neutron star is not expected 
to be perfectly Planckian, however (van Paradijs 1982; 
London et al. 1984, 1986; see also Titarchuk 1994; Madej 
1997, and references therein). This is mainly due to the 
effects of electron scattering in the neutron star atmo- 
sphere, deforming the original X-ray spectrum. This re- 
sults in a systematic difference between the effective tem- 
perature (as would be measured on Earth), T e ff iOCJ , and 
the temperature as obtained from the spectral fits, Tbb 
(also referred to as 'colour' temperature, see e.g. Lewin 
et al. 1993). In general, the deviations from a Planckian 
distribution will depend on several parameters, such as 
temperature, elemental abundance, neutron star mass and 
radius. The hardening factor, Tbb/T e ff,oo ; has been deter- 
mined through numerical calculations by various people 
and its value is typically around 1.7. When the burst lu- 
minosity approaches the Eddington limit the deviations 
from a black-body become larger, and so does the spec- 
tral hardening (Tbb/T c ff i00 ^2, Babul & Paczyhski 1987; 
Titarchuk 1988). During extreme radius expansion phases, 
however, this trend may break down and Tbb/T e ff j0 o<l 
(Titarchuk 1994). Attempts have been made to determine 
the spectral hardening from the observed cooling tracks, 
but conclusions are still rather uncertain (e.g. Penninx et 
al. 1989; Smale 2001). As a result, the interpretation of X- 
ray bursts spectra has remained uncertain and constraints 
on the mass-radius relationship for neutron stars elusive. 

Type I X-ray burst theory predicts three different 
regimes in mass accretion rate (M) for unstable burn- 
ing (Fujimoto et al. 1981, Fushiki & Lamb 1987; see also 
Bildsten 1998, 2000, Schatz et al. 1999, and references 



therein; note that values of critical M depend on metal- 
licity, and on assumed core temperature and mass of the 
neutron star): 

1) low accretion rates; 10~ 14 M yr" 1 <; M <5 2 x 
10~ 10 Mq yr _1 : mixed H/He burning triggered by 
thermally unstable H ignition 

2) intermediate accretion rates; 2 x 10~ 10 Mq yr~ 4 ^ 
M <; 4-11 x 10~ 10 Mq yr" 1 : pure He shell ignition after 
steady H burning 

3) high accretion rates; 4-11 x 10~ 10 Mq yr _1 <J M <, 
2 x 10~ 8 Mq yr _1 : mixed H/He burning triggered by 
thermally unstable He ignition 

H and He are burning stably in a mixed H/He environ- 
ment for very low and very high values of M, i.e. M be- 
low ~10~ 14 M yr" 1 and above -2 x 10~ 8 M yr" 1 (close 
to the critical Eddington M). During pure helium flashes 
the fuel is burned rapidly, and such bursts therefore last 
only 5-10 s. This gives rise to a large energy release in 
a short time, which causes the bursts often to reach the 
Eddington limit, leading to photosphcric radius expan- 
sion. Bursts with unstable mixed H/He burning release 
their energies on a longer, 10-100 s, timescale, due to the 
long series of (3 decays in the rp-process (see e.g. Bildsten 
1998, 2000). 

The Z sources (Hasinger & van der Klis 1989) are a 
group of sources inferred to persistently accrete near the 
Eddington limit. In a colour-colour diagram they trace 
out a Z-like shape, with the three limbs of the Z (his- 
torically) referred to as the horizontal branch (HB), nor- 
mal branch (NB) and flaring branch (FB), from top to 
bottom. M is inferred to increase from sub-Eddington at 
the HB, near-Eddington at the NB to super-Eddington at 
the FB (e.g. Hasinger 1987; Lamb 1989; Hasinger et al. 
1990). According to the burning regimes outlined above 
these sources should exhibit long (^ 10-100 s) type I X- 
ray bursts, at least on the HB and NB. However, of 
the Z sources, only GX 17+2 and CygX-2 show (infre- 
quent) bursts (Kahn & Grindlay 1984; Tawaraet al. 1984c; 
Sztajno et al. 1986; Kuulkers et al. 1995, 1997; Wijnands 
et al. 1997; Smale 1998), indicating that most of the ma- 
terial is burning stably. This is in contrast to the above 
theoretical expectations for high M. Moreover, the bursts 
in CygX-2 are short (—5s), whereas GX17+2 shows both 
short (~10s) and long (^100s) bursts. One of the bursts 
of CygX-2 showed a radius expansion phase (Smale 1998); 
this burst clearly bears all the characteristics of a He flash 
(regime 2), whereas the neutron star is inferred to ac- 
crete at near-Eddington rates. The short duration bursts 
in GX 17+2 also hint to a He flash origin, whereas the 
long duration bursts hint to unstable mixed H/He burn- 
ing (regime 3; see also van Paradijs et al. 1988). A study 
of X-ray bursts from CygX-2 and GX17+2 observed by 
EXOSAT showed no correlation of the burst properties 
with position in the Z, although the number of bursts 
observed from GX17+2 was small (Kuulkers et al. 1995, 
1997). 
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The Proportional Counter Array (PCA) onboard the 
Rossi X-ray Timing Explorer (RXTE) combines high 
throughput using a large collecting area (maximum of 
~6500cm 2 ) with the ability to label photons down to a 
time resolution of /is. This is ideal to study short events 
like X-ray bursts, especially during the start of the burst. 
Such studies may provide more insight in the properties 
of X-ray bursts which occur at these extreme mass accre- 
tion rates. Analysis of X-ray bursts in persistent sources 
at (relatively) high inferred mass accretion rates (typi- 
cally £0.2 M Ed d) observed with the RXTE/PCA were pre- 
sented for one burst seen with CygX-2 (Smale 1998) and 
one seen in GX3+1 (Kuulkers & van der Klis 2000). In 
this paper we present the first account of ten X-ray bursts 
from GX17+2 observed by the RXTE/PCA during the 
period 1996-2000. For a description of the correlated X- 
ray timing and spectral properties of GX 17+2 using the 
same data set we refer to Homan et al. (2001). 

2. Observations and Analysis 

The PCA (2-60 keV; Bradt et al. 1993) onboard RXTE 
observed GX 17+2 various times during the mission. Up 
to now a total of 657ksec of useful data has been obtained. 
A log of these observations is given in Table 1. During 
the observations in 1996-1998 all five proportional counter 
units (PCUs) were active, whereas in 1999 and 2000 only 
three units were active. The high voltage settings of the 
PCUs have been altered three times during the mission 
(so-called gain changes), which modified the response of 
the detectors. These changes therefore mark four periods, 
called gain epochs. The observations were done in two 
standard modes: one with relatively high spectral resolu- 
tion (129 energy channels covering the whole PCA energy 
band) every 16 sec, the Standard 2 mode, the other hav- 
ing no spectral information providing the intensity in the 
whole PCA energy band at a moderate time resolution of 
0.125 s, the Standard 1 mode. Additionally data were 
recorded in various high time resolution (<2ms) modes 
that ran in parallel to the Standard modes, and that 
recorded photons within a specific energy band with ei- 
ther low spectral resolution (B-modes or E-modes) or no 
spectral resolution (SB-modes). The B-, E- and SB-modes 
used here combined the information from all layers of all 
active PCUs together. During the 1996 observations a B- 
and E-mode were available, giving 16 and 64 energy bands, 
covering channels 0-49 and 50-249, at 2 ms and 125 /is, re- 
spectively. For most of the observations in 1997 four SB 
modes covering the total PCA energy range were available. 
In the 1998, 1999, and 2000 observations two SB-modes 
covering channels 0-13 and 14-17 at 125/xs, and an E- 
mode giving 64 energy bands covering channels 18-249 at 
16 /is time resolution were available. 

For the spectral analysis of the persistent emission 
we used the Standard 2 data. We accumulated data 
stretches of 96 s just before the burst, combining the PCUs 
which were operating at that time. In order to study the 
spectral properties of the bursts we used two approaches. 




5000 10 4 1.5x10 4 2x10 4 

Time (sec) 

Fig. 1. Standard 1 light curve of the GX17+2 observa- 
tions during October 10, 1999, at a time resolution of 5 s. 
Time zero corresponds to 05:39:27 (UTC). No corrections 
for background and dead-time have been applied. Clearly, 
the source varies on the same time scale (and faster) as the 
burst which started at T = 12680 s. The source was in the 
FB and the lower part of the NB during the observations. 

Spectra during the bursts were determined every 0.25 sec 
for the first ~20s of the burst. For the short bursts this 
means the whole duration of the burst. For the long bursts 
we also used the Standard 2 data to create spectra at 
16 s intervals, for evaluating the remainder of these bursts. 
Since no high time resolution spectral data were available 
during the 1997 observations (only 4 SB-modes), only the 
Standard 2 data were used to study the spectral prop- 
erties of the burst from this observation. All spectra were 
corrected for background and dead-time using the proce- 
dures supplied by the RXTE Guest Observer Facility^. 
A systematic uncertainty of 1% in the count rate spectra 
was taken into account. For our spectral fits we confined 
ourselves to the energy range of 3-20 keV, which is best 
calibrated. The hydrogen column density, Nh, towards 
GX17+2 was fixed to that found by the Einstein SSS 
and MPC measurements (2xl0 22 atoms cm -2 , Christian 
& Swank 1997; see also Di Salvo et al. 2000). In all cases 
we included a Gaussian line (see Di Salvo et al. 2000) fixed 
at 6.7keV, with a fixed line width of 0.1 keV. One sigma 
confidence errors were determined using A% 2 = 1. 

Large amplitude, high coherence brightness oscilla- 
tions have been observed during various type I X-ray 
bursts in other low-mass X-ray binaries (LMXBs; see e.g. 
Strohmayer 1998, 2001). We searched all the bursts from 
GX17+2 for such oscillations. Using the high time res- 
olution modes we performed Fast Fourier Transforms to 

1 http: / /heasarc. gsfc.nasa.gov /docs / xte /recipes/cook_book.html. 



4 E. Kuulkers et al.: X-ray bursts in GX 17+2 

Table 1. RXTE observation log of GX17+2" 



Year Start (UT) End (UT) #°T 

(ksec) bursts 



1996 


feb 


07 


13:27 


feb 


09 


00:02 


58 


1 


1997 


feb 


02 


19:13 


feb 


27 


03:34 


59 


1 


1997 


apr 


01 


19:13 


apr 


04 


23:26 


35 





1997 


Jul 


27 


02:13 


Jul 


28 


00:33 


43 





1998 


aug 


07 


06:40 


aug 


08 


23:40 


71 


1 


1998 


nov 


IS 


06:42 


nov 


20 


13:31 


86 


2 


1999 


oct 


03 


02:43 


oct 


12 


07:05 


298 


5 


2000 


mar 


31 


12:15 


mar 


31 


16:31 


7 






a Between Feb 2-27, 1997 observations every 3-6 days. 
6 Total effective exposure time. 



produce power spectra with a Nyquist frequency of 2048 
Hz. This was done in the total energy band (2-60 keV) and 
in a high energy band ranging from ~8keV to ~20keV. 
For the 1996 observations, however, the high time reso- 
lution mode only covered the 13.5-60 keV range; the 2- 
60keV power spectra had a Nyquist frequency of 256 Hz. 
In all cases, our searches were carried out on power spec- 
tra of 0.25 s or 2 s in duration. To increase the sensitivity 
for cases where burst oscillations are only present for a 
period of time comparable or shorter than the length of 
the power spectrum, we 'oversampled' the data by factors 
of 2 and 8, respectively, by taking the start time of the 
next data segment to be 0.125 s and 0.25 s later than that 
of the previous one, instead of 0.25 s and 2 s (i.e. we use 
overlapping data segments) . 

3. Results 

3.1. Temporal behaviour 

Inspection of all the Standard 1 light curves of GX 17+2 
at 1 s time resolution yielded a total of fourteen burst- like 
events (designated bl, b2, blO, fl, f4; the reason 
for different designations is the nature of the events, see 
below). Seven events were rather short, ~10s (bl, b3, b5, 
fl-f4), whereas the other events were long, ~360-1600s 
(b2, b4, b6-bl0). Table 2 gives the start times of these 
events, together with the gain epoch in which the obser- 
vations were done and the number of active PCUs. In that 
table we also give the duration of the events as estimated 
by eye using the light curves in the total PCA energy 
band and/or the hardness (ratio of count rates in two en- 
ergy bands) curves, the rise time of the event as deter- 
mined in two ways (see below), and the exponential decay 
time as determined from a fit to the decay part of the 
Standard 1 light curves, together with the goodness of 
this fit expressed in terms of reduced x 2 , X 2 c d- We note 
that the estimates of the durations are rather arbitrary, 
but robust methods such as the so-called T90 measure- 
mentsf] as used for 7-ray bursts (see Koshut et al. 1996, 
and references therein) can not be applied, especially to 

2 T90 is defined as the time it takes to observe 90% of the 
total background-subtracted counts in an event, starting and 



the longer events, since the persistent emission varies on 
the same time scale (or even fatser) as the event itself, see 
e.g. Fig. [I]. The rise time of the events was determined as 
follows. We constructed light curves in the full PCA en- 
ergy band with a time resolution of 1/32 sec. We defined 
t r i SO as the time for the event to increase from 25% to 90% 
of the net peak event rate (see e.g. Muno et al. 2000; van 
Straaten et al. 2001). A detailed look at the light curves 
revealed that the onset of events bl-b4, b8, blO, and f2 
consisted of a rapid rise phase and a subsequent slower 
rise to maximum. In the events b5-b7, b9, fl, and f3 the 
count rate rapidly increased to maximum, with no sub- 
sequent slower rise. We therefore also fitted the pre-event 
phase, fast rise phase (and slow rise phase) with a constant 
level and one (or two) linear functions, respectively, using 
the Standard 1 light curves. We find that the total du- 
ration of the fast rise phase, tf r , is between 0.1 and 0.6 sec 
for events bl-blO and 0.4-1 s for fl-f4 (Table 2). In the 
fits to the decay portion of the light curves for the short 
events we used the Standard 1 light curves at 0.125 sec 
time resolution, while for the longer events we rcbinncd 
these light curves to a time resolution of 5 sec. For some 
of the events an exponential does not describe the decay 
very well; this is probably due to the short term variations 
in the persistent emission. For these events we fitted only 
the initial decay (first few seconds for f2 and f3, and first 
few 100s for burst b9 and blO). 

In Fig. H we show the light curves of the four events 
fl-f4, at low (<J7keV) and high (^7keV) energies, with 
the corresponding hardness (ratio of the count rates in the 
high and low energy band) curves, all at a time resolution 
of 0.125 s. Although they have a fast rise and a longer 
decay (see also Table 2), they show small or no variations 
in hardness. Time-resolved spectral analysis (like done in 
Sect. 3.2.3) confirmed this; no clear cooling during the 
decay can be discerned. We can therefore not classify these 
events as type I bursts. Since all of the four events occurred 
in the FB, we conclude that they must be flares, of which 
the light curves happen to resemble those of X-ray bursts 
(such as bl, b3 and b5). We will not discuss these four 



ending when 5% and 95% of the total background-subtracted 
counts have been observed. 



E. Kuulkers et al.: X-ray bursts in GX 17+2 
Table 2. Bursts and burst-like events (flares) in GX17+2 



event" 




Start Time (UT) 


1)1 




feb 


08 


16:17:12 


b2 


1997 


fel) 


08 


02:36:34 


b3 


1998 


aug 


07 


13:15:50 


b4 


1998 


nov 


18 


08:51:26 


b5 


1998 


nov 


18 


14:37:30 


b6 


1999 


oct 


03 


15:36:32 


b7 


1999 


oct 


05 


23:41:43 


b8 


1999 


oct 


06 


11:10:33 


b9 


1999 


oct 


09 


12:34:24 


blO 


1999 


oct 


10 


09:10:47 


ft 


1996 


fel) 


07 


03:39:11 


f2 


1998 


nov 


19 


14:38:24 


f3 


1998 


nov 


20 


00:44:43 


f4 


1999 


oct 


11 


08:55:30 



E" 



PCUs c 
5 
5 
5 
5 
5 

(0,2,3) 
(0,2,4) 
(0,2,3) 
(0,2,3) 
(0,2,3) 

5 
5 
5 

3 (0,2,3) 



tdu 



t 9 



(rcd/dof' 1 branch 1 



10 

>360 
10 

1000 
10 

1600 
500 
500 
500 
700 

10 
10 
10 
10 



1.22 
1.19 
0.53 
1.34 
0.41 
0.41 
0.56 
1.66 
0.13 
0.72 

1.03 
0.41 
0.75 
5.0 



0.35 ± 0.05 
0.34 ± 0.08 
0.27 ± 0.09 
0.61 ± 0.04 
0.54 ± 0.04 
0.19 ± 0.04 
0.30 ± 0.03 
0.13 ± 0.02 
0.16 ± 0.02 
0.20 ± 0.04 

1.12 ± 0.07 
0.33 ± 0.04 
0.45 ± 0.06 
1.06 ± 0.20 



1.83 ± 0.08 
248 ±J 

2.55 ± 0.24 
197 ± 2 

2.06 ±0.13 
274 ± 3 

77.3 ± 1.2 
70.2 ± 1.4 

76.4 ± 1.5 
115 ± 3 

2.98 ± 0.49 
1.72 ± 0.35 
1.85 ± 0.25 
2.18 ± 0.22 



1.1/130 
2.5/49 
1.0/131 
3.7/147 
1.1/85 
2.0/242 
1.9/104 
1.2/57 
3.1/66 
3.7/57 

1.3/88 
1.1/28 
1.1/55 
1.1/148 



mNB 
mNB 

1NB 

SV 
mNB 

1HB 
uNB 

1HB 
uNB 

1NB 

mFB 
uFB 
uFB 
1FB 



a Designation used in text. 

6 RXTE gain epoch of the observation. 

c Numbers of active PCUs; if <5 the PCU unit numbers are given between brackets. 
d Approximate duration in sec. 

e Time (sec) for the count rate to rise from 25% to 90% of the net-peak rate. 
' Duration of fast rise phase in sec, see text. 
9 Exponential decay time in sec. 

h Goodness of the exponential fit to the decay part of the light curve. 

1 Source position in the Z just before the event (see also Fig. |l6|); u, m, and 1 stand for upper, middle and lower, 
respectively. SV means 'soft' vertex, i.e. the NB/FB transition region. 



f 1 



«8 %d 





f3 







10 10 10 

Time after start of burst (sec) 

Fig. 2. The X-ray flare light curves at low (upper panels) and high (middle panels) energies and the corresponding 
hardness curves (lower panels), all at a time resolution of 0.125 sec, for the four flares (fl— f4; see text). Time zero is 
the flare start time as given in Table 2. The low and high energy ranges are 1.5-7.2 keV and 7.2-19.7keV, respectively, 
for fl, 1.9-6.2keV and 6.2~19.6keV, respectively, for f2 and f3, and 2.1-7.1keV and 7.1-19.9keV, respectively, for 
4. Hardness is defined as the count rate ratio of the high to low energy band. No corrections for background and 
dead-time have been applied. 
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CO 

d 



b3 







o 



Time after start of burst (sec) 



Fig. 3. Same as Fig. |[ but for the three short X-ray bursts bl, b3, and b5. The low and high energy ranges are 
l-7.2keV and 7.2-19.7keV, respectively, for bl, whereas they are 1.9-6.2keV and 6.2-19.6keV, respectively, for b3 
and b5. 



events further in the paper and denote the remainder of 
the ten events as bursts, since we will show below that 
they are genuine type I X-ray bursts. 

In Fig. |^ we show the light curves of the three short 
bursts bl, b3 and b5, at low and high energies, with the 
corresponding hardness curves, all at a time resolution of 
0.125 s. All three bursts show a fast rise (typically less than 
0.5 s) and an exponential decay with a decay time of ~2s 
(see Table 2). During the rise the emission hardens; as 
the bursts decay, the emission becomes softer. The main 
difference between the three bursts is the fact that the 
peak of burst b3 is about 25% lower than the other two 
bursts; it looks like a 'failed' burst. It also has two peaks, 
as if some new unstable burning occurred, ~5 s after the 
start of the burst. 

In Figs. [| and || we show the light curves of the long 
bursts b2, b4, and b6-bl0, at low and high energies, with 
the corresponding hardness curves, all at a time resolu- 
tion of 2 s. In Figs. || and ^ we focus on the start of these 
bursts, all at a time resolution of 0.125 s. Again the rise 
times are very short (also typically less than 0.5 s), but the 
decay times are much longer, with decay times in the range 
~70-280s (see Table 2). Apart from the fact that the hard 
burst emission decays faster than the soft burst emission 
(i.e. spectral softening), there are more pronounced differ- 
ences between the light curves in the two energy bands. 
In Figs. f| and || one sees that all the low energy light 
curves show a kind of spike at the start of the decay. These 
spikes last for a few seconds in most cases; however, for 
burst b6 it seems to last ~15s (with an exponential de- 
cay time of 9±2s). At high energies no such spikes occur 



(except for burst blO); instead the bursts have more flat- 
topped peaks, with durations ranging from tens of seconds 
to ~200s. Burst b6 is the nicest example of this. Zooming 
in on the start of these long bursts, it becomes clear that 
the rise is somewhat slower at high energies with respect 
to low energies (causing the hardening of the emission dur- 
ing the early phase of the burst). Also, in bursts b4 and 
b6-b9 very short (<0.5s) spikes occur during the rise in 
the high energy light curve, which again are especially ev- 
ident in burst b6 (two spikes!). This causes the hardness 
to drop on the same time scale in some of these bursts. 
Later we will show that this corresponds to fast radius 
expansion/contraction episodes. 

3.2. Spectral behaviour 

3.2.1. Persistent emission before the bursts 

The persistent emission just before the bursts b6 and b7 
can be satisfactorily (x;? ed =0.7 and 1.2, respectively, for 35 
degrees of freedom, dof) described by an absorbed cut-off 
power-law component plus a Gaussian line. For the per- 
sistent emission before the remainder of the bursts this is 
not a satisfactory model (x^ed ranges from 1.4 with 35 dof 
for b8 to 5.5 with 42 dof for b3). An additional compo- 
nent is necessary to improve the fits. We used the F-test 
to calculate whether the additional component was indeed 
significant. For the additional component we chose a black 
body, as is generally used when modeling the X-ray spec- 
tra of bright LMXBs (e.g. White et al. 1986; Christian 
& Swank 1997; Church & Balucihska-Church 2001). Note 
that more complicated models are necessary to describe 
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Fig. 4. Same as Fig. [| but at a time resolution of 2 s, for the long bursts b2, b4 and b6. Note that b2 was interrupted 
by a South Atlantic Anomaly (SAA) passage as can be seen by the sudden decrease in count rate. The low and high 
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and 2.1-7.1 keV and 7.1-19.9 keV, respectively, for b6. 



Table 3. Persistent emission spectral parameters 



# 


F b 

- 1 - pcrs 




kT bb 


Rbb,10 C 




E cu t 


norm. e 


gnorm/ 


Xred / 


bb 9 








(2-20 keV) 




(keV) 


(km) 




(keV) 






dof 








bl 


2.0±0.3 


1 


15±0.06 


14.1±1.5 


1.0+0.1 


4.5+0.2 


3.8±0.4 


0.013+0.002 


0.89/47 


16 


5x10" 


10 


b2 


2.2+0.5 


1 


10±0.03 


22.1±2.2 


0.4+0.2 


3.8±0.2 


1.5±0.4 


0.008±0.002 


0.56/40 


29 


4x10" 


18 


b3 


1.9±0.3 


1 


13±0.03 


19.2±1.5 


0.9±0.1 


4.2+0.2 


2.6+0.4 


0.013+0.002 


0.58/40 


29 


1x10" 


18 


b4 


2.0±0.3 


1 


10±0.04 


20.3±1.8 


0.7±0.1 


4.0±0.2 


2.3+0.4 


0.014±0.002 


0.85/40 


28 


6x10" 


15 


b5 


2.3±0.3 


1 


16±0.05 


17.5±1.7 


0.9+0.1 


4.5±0.2 


3.2+0.4 


0.010±0.002 


0.59/40 


23 


5x10" 


15 


b6 


2.5±0.1 








1.03±0.03 


5.1±0.1 


4.8±0.1 


0.012±0.002 


0.71/35 


<r 


8 




b7 


2.5±0.1 








1.17±0.03 


4.6±0.1 


7.0±0.2 


0.015±0.002 


1.22/35 


<2 i 


30 




b8 


2.4±0.5 


1 


08+0.13 


11.0±4.3 


1.0+0.1 


5.3+0.3 


3.8+0.5 


0.011±0.002 


1.13/33 


6 


1 




b9 


2.4+0.4 


1 


26+0.08 


11.6±1.5 


1.0+0.1 


4.9+0.3 


4.3+0.4 


0.010±0.003 


0.87/33 


13 


8x10 


-5 


blO 


1.9±0.3 


1 


18±0.05 


15.1±1.4 


1.1+0.1 


4.5+0.3 


3.7±0.5 


0.011±0.002 


1.16/33 


22 


5x10- 


-7 



a Nh was fixed to 2xl0 : 
6 



atoms cm 



see text. 

Unabsorbed persistent flux in 10 -8 ergs -1 cm" 2 . 
c Apparent black-body radius at lOkpc. 
d Power-law photon index. 

e Power-law normalization in photons keV~ cm - s _1 at 1 keV. 

* Normalization in photons cm -2 s _1 in the line modeled by a Gaussian fixed at 6.7 keV and width (a) of 0.1 keV. 
9 Black-body emission contribution to the total persistent flux in % (2-20 keV). 
h Probability (in%) that inclusion of a black-body component is not significant. 
1 95% confidence upper limit, using kTbb = l-14 keV, see text. 



the spectra of GX 17+2 with better energy resolution and 
broader energy range (e.g. Di Salvo et al. 2000). However, 
our approach is adequate for this paper. In the rightmost 
column of Table 3 we give the results of our F-tests, i.e. 
the probability that the addition of the black-body com- 
ponent in the spectral fits is not significant. We conclude 
that except for bursts b6 and b7, an additional component 



is warranted on the <;99% level. In Table 3 we present the 
results of the spectral fits to the persistent emission. Some 
of the x 2 od values are rather small (bursts b2, b3, b5), per- 
haps due to a slight overestimate of the systematic error 
added (1%). The unabsorbed persistent flux just before 
the bursts during the various observations varied between 
1.9-2. 5xl0 -8 ergs -1 cm -2 (2-20 keV). For bursts b6 and 
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Fig. 5. Same as Fig. ||, for the last four long bursts observed in 1999 Oct (b7-bl0). The low and high energy ranges 
are 1.9-7.1keV and 7.1-19. 9keV, respectively. 
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Fig. 6. Same as Fig. |[ but 5 s before and 10 s after the 
start of the bursts b2, b4 and b6, at a time resolution of 
0.125s. 



b7 we determined single parameter 95% confidence up- 
per limits (using Ax 2 =2.71) on the black-body contribu- 
tion, by including a black-body component in the spec- 
tral fits, and fixing the temperature to its mean value de- 
rived for the persistent spectra of the other bursts, i.e. 
kTbb=l-14keV. The black-body component contribution 



to the persistent emission varied between less than 2% 
(burst b7) up to 29% (bursts b2, b3) in the 2-20 keV band 
(Table 3). 



3.2.2. Persistent emission during the bursts? 
3.2.2.1 Previous EXOSAT results 

Usually it is assumed that the persistent emission is not 
influenced by the burst and that one can, therefore, study 
the burst by subtracting the persistent emission from the 
total source emission. This is referred to as the 'stan- 
dard' burst spectral analysis (see e.g. Sztajno et al. 1986). 
However, if the neutron star photosphere contributes sig- 
nificantly to the persistent emission, this approach is not 
correct, if the burst emission originates from the same 
region (van Paradijs & Lewin 1986). In this case the spec- 
tral fits to the net burst spectra yield a systematically 
larger black-body temperature, Tbb, and smaller appar- 
ent black-body radius, Rbb, especially near the end of 
the burst when the net burst flux is low. In fact, in this 
case of incorrect subtraction of the persistent emission, 
the net burst spectrum is not a black-body. Van Paradijs 
& Lewin (1986) proposed to fit the total source spectrum 
with a two-component model, a black-body component 
and a non-black-body component. During the burst the 
non black-body component is fixed to what is found in 
the persistent emission, and the black-body component 
is left free. The black-body component should include all 
emission from the neutron star photosphere. The underly- 
ing idea is that the non black-body component arises from 
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Fig. 8. Leftmost panel: Two-component spectral fit results for the total burst emission of burst b4 plotted on a 
logarithmic time scale. The filled circles and open squares represent the fit results of the 0.25 s and 16 s spectra, 
respectively. The data have been logarithmically rebinned in time for clarity. The values for the persistent black-body 
component have been indicated by a dotted line. Right panels: Two-component spectral fit results for the 0.25 s spectra 
of the total burst emission of burst bl, b4 and b6, see text. The values for the persistent black-body component have 
been indicated by a dotted line for burst bl and b4. Both panels: from top to bottom: bolomctric black-body flux, Fbt> ? 
in 10~ 8 ergs -1 cm -2 , black-body temperature, kTbb, apparent black-body radius, Rbb.io, at lOkpc, and goodness of 
fit expressed in Xred- For bursts bl, b4 and b6 the number of dof is 22, 18 and 16, respectively, for the 0.25s spectral 
fits. For burst b4 the number of dof is 44 for the 16 s spectral fits. Note the difference in scales of Xrcd m the leftmost 
panel with respect to the other panels. 
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Fig. 9. a: At the top the first 16 s spectrum observed after the start of burst b2 is displayed. A two-component fit is 
shown, i.e. a single black-body and cut-off power-law plus Gaussian line (subjected to interstellar absorption). The 
parameters of the cut-off power-law component plus Gaussian line are fixed to the values derived for the persistent 
emission. At the bottom the residuals after subtracting the best model from the observed spectrum is displayed. The fit 
is clearly bad (x re( j/dof=24.2/44). b: At the top the same first 16s spectrum at the beginning of burst b2 is displayed. 
Now a three-component fit is shown, i.e. two black-body components and one cut-off power-law component (subjected 
to interstellar absorption). The cut-off power-law and one black-body component parameters are fixed to the values 
derived for the persistent emission. At the bottom the residuals after subtracting the best model from the observed 
spectrum is displayed. The fit has clearly improved (x^ ed /dof=1.8/44). 



the accretion process, and is not influenced by the X-ray 
burst. 

GX 17+2 is a bright X-ray source, and presumably 
the neutron star contributes significantly to the persistent 
emission (Van Paradijs & Lewin 1986; Sztajno et al. 1986). 
Since the source is bright, compared to most other burst 
sources the net flux at the peak of the burst relative to the 
persistent flux is rather low. Sztajno et al. (1986) found 
that the black-body component contributed ~40% to the 
persistent emission just before the two bursts observed 
by EXOSAT. Using the 'standard' burst spectral analysis 
Sztajno et al. (1986) found relatively high black-body tem- 
peratures (kTbb— 2-3 keV) and relatively small apparent 
black body radii at a distance of lOkpc (Rbb.io— 3-5 km, 
for isotropic emission) for the short (~10s) burst. For the 
long (>5min) burst, kTbb only showed a small change 
from ~2.1keV at the peak of the burst to ~1.7keV near 
the end of the burst, with Rbb,io decreasing from ~7km 
to ~4km. Their fits were, however, satisfactory, with x re d 
of 0.6-1.2. Using the two-component model, Sztajno et 
al. (1986) found that during the short burst the values 
for Tbb were in the range normally seen in type I X-ray 
bursts, and that the systematic decrease in Rbb had disap- 
peared. The Xred va hies for the two-component model fits 
ranged between 0.6 and 1.3, except for one fit, for which 
Xred = l-6, i-e. slightly worse than the 'standard' spectral 
analysis. They attributed this to the smaller relative er- 
ror in the total burst data than in the net-burst data, for 
which the persistent emission was subtracted. Note that 
in this case Sztajno et al. (1986) found a slight increase 



in Rbb, io, apparently anti-correlated with Tbb 7 which they 
argue was due to the non-Planckian shape of the spectrum 
of a hot neutron star (van Paradijs 1982; see Titarchuk 
1994; Madej 1997, and references therein; see, however, 
Sect. 5.2). 

3.2.2.2 Our RXTE results 

Guided by the results of Sztajno et al. (1986) discussed 
in Sect. 3.2.2.1, we first fitted the total burst data using 
a black body and a cut-off power law plus a Gaussian 
line. The parameters of the absorbed cut-off power law 
and Gaussian line were fixed to the values found for the 
persistent emission before the burst (Table 3). Using this 
model we obtained good fits to the 16 s spectra of bursts 
b6 and b7, for which the persistent emission did not con- 
tain a significant black-body contribution (x re d 01 —1 for 
37 dof). However, the 16 s spectral fits were bad when- 
ever the persistent emission spectra contained a black- 
body component; the fits became worse as the persis- 
tent black-body contribution became stronger (for the 
first -100s of the burst: xLi/dof- 1 - 5 - 3 / 37 [ burst b8 L 
X rcd /dof~2. 5-8/37 [burst b9], x? d/ dof - 5 - 10 / 37 I burst 
b!0], x?ed/ dof - 2 °- 24 / 44 I burst b4 L X? cd / dof - 22 - 25 / 44 
[burst b2] ) . The worst x TO d, occurred near the peak of the 
bursts. For instance, Fig. || shows the best-fit parameters 
and x?od f° r bur st b4, for which the persistent emission 
had a black-body contribution of ~28% (2-20 keV). An 
example of a burst spectrum and the result of the two- 
component fit is shown in Fig. ^a for burst b2. The Xred 
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Fig. 10. a: Bolometric peak black-body flux (see text), 
Fbb.peak, in units of 10~ 8 ergs -1 cm -2 versus the bolomet- 
ric persistent black-body flux, Fbb.pers- Radius expansion 
(RE) bursts are denoted with an open square. The dotted 
line represents a fit to the data points (excluding burst b2 
and b3) to the function Fbb.pcak = Fbb.pcrs + C. The con- 
stant C corresponds to (1.31±0.03)xl0 _8 ergs _1 cm" 2 . 
b: Bolometric peak black body flux from spectral fits 
after subtraction of the persistent emission (see text), 
F b b,not_pcak, in units of 1CT 8 ergs -1 cm -2 versus F b b, P crs- 
Radius expansion (RE) bursts are denoted with an open 
square. The dotted line represents the constant C. 



with 44 dof for burst b4 decreases from ~24 at the be- 
ginning of the burst to ~2 near the end of the burst. A 
similar trend is seen for the 0.25 s burst spectral fits: Xred 
of ~2 for 18 dof during the first ~10s of burst b4 and an 
X 2 ed of ~1 for 16 dof during the first ~10s of burst b6 
(see Fig. ||). The x 2 ed values for 0.25 s spectral fits were 
much lower than those for the 16 s spectral fits, due to the 
much lower signal to noise and lower spectral resolution 
of the 0.25 s spectra. Note that the high values of Xred 
for the 16 s spectra are not due to fast spectral variations 
within the time the spectra are accumulated. This only 
applies to spectra which include the first few seconds of 
the bursts, during the short radius expansion and initial 
contraction phase (see Sect. 3.2.4). Halfway the decline of 
burst b4, Rbb seems to increase again; this is also seen in 
the other long bursts (except burst b6), and in the short 
bursts (e.g. burst bl, see Fig. ||). This is similar to that 
seen by Sztajno et al. (1986), but their fits seemed better 
(see Sect. 3.2.2). This is (probably) due to the much lower 
signal to noise ratio of the EXOSAT/ME data. 

The mean bolometric black-body flux during the 
first 10 s after the peak, Fbb.pcak, of bursts b4 is 



higher than that of burst b6. The difference is 
0.63(±0.12)xl0 _8 ergs _1 cm~ 2 (for the uncertainty in 
Fbb.pcak we used the rms variation in the bolometric 
black-body fluxes). This is consistent with the bolo- 
metric black body contribution in the persistent emis- 
sion, Fbb.persi before burst b4 (indicated by the dot- 
ted lines in the top panels for burst b4 in Fig. ||): 
Fbb,pcrs=0.65(±0.05)xl0~ 8 ergs~ 1 cm -2 . We measured 
Fbb.pcak and Fbb,pcrs for all the bursts to see if this ef- 
fect is seen in other bursts as well. We used the highest 
observed black body flux, except for burst b2 where we 
used the first 16 s measurement after the start of the burst. 
The result is displayed in Fig. [io|a. Indeed, Fbb.pcak differs 
between the bursts and is clearly correlated with Fbb.pers 
(except for burst b3). Such a relation can be expected 
if Fbb.peak is close to a certain upper limit (presumably 
the Eddington limit) and includes a (certain fraction) of 
the persistent emission. We can reasonably fit the data 
points to the function Fbb.peak = Fbb.pers + C (excluding 
bursts b2 and b3; burst b2 did not cover the first 10 s after 
the peak of the burst, and burst b3 was weaker than the 
other bursts). The resulting x 2 od =2.0 for 7 dof. We find 
C = (1.31 ± 0.03) xl0~ 8 ergs" 1 cm" 2 ; this is close to the 
values of Fbb.peak for bursts b6 and b7 which had no sig- 
nificant black-body contribution in the preceding persis- 
tent emission. In Fig. [l(] we have also indicated whether 
a burst was a radius expansion event (see Sect. 3.2.4). 
It can be seen in Fig. |lO| a that the bursts with no ra- 
dius expansion/contraction phase (except burst b3) have 
somewhat higher values of Fbb.peak than the radius ex- 
pansion bursts (except burst b4). Theoretically, during 
the radius expansion/contraction phases the flux should 
equal the Eddington limit. Therefore, the (peak) fluxes 
observed for different radius expansion bursts should have 
similar values, whereas for bursts with no radius expan- 
sion/contraction phase they should be either similar or 
smaller. This is not the case. 

To summarize, we find that the method proposed by 
Van Paradijs & Lewin (1986) and applied by Sztajno et al. 
(1986) to the bursts observed by EXOSAT does not work 
for our bursts observed with the RXTE/PCA. The two- 
component spectral fits during the bursts give bad results 
whenever there is a black-body contribution to the per- 
burst persistent emission. Moreover, the total burst peak 
fluxes, Fbb.peak are different from burst to burst and de- 
pend on the amount of the black-body contribution to the 
persistent emission, Fbb.pers- This is not what one would 
expect if during some of the bursts a limit is reached (pre- 
sumably the Eddington limit). We conclude that the per- 
sistent black-body emission does not disappear during the 
burst. An example is displayed in Fig. ^b. We show again 
the first 16 s spectrum during burst b2 plus now the results 
of a three-component fit, i.e. two black bodies and a cut- 
off power-law plus Gaussian line (subjected to interstellar 
absorption). The parameters of the cut-off power law plus 
Gaussian line and one black-body component have been 
fixed to those found for the persistent emission (so the fit 
does have the same number of dof as the one displayed 
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Fig. 11. Spectral fit results for the net burst emission 
of burst bl, b3 and b5; from top to bottom: bolometric 
black-body flux, Fbb,not, in 10~ 8 ergs -1 cm~ 2 , black-body 
temperature, kTbb, apparent black-body radius, Rbb.io, at 
lOkpc, and goodness of fit expressed in reduced x 2 ■ The 
number of dof is 22 for burst bl and 18 for bursts b3 and 
b5. 



in Fig. ^p,). Clearly, the fit has improved considerably. 
This means that the burst emission is decoupled from the 
persistent emission, and therefore the 'standard' spectral 
analysis should work in this case. This is the subject of 
the next subsection. 

3.2.3. 'Standard' spectral analysis: net burst emission 

When modeling the net burst spectra (i.e. spectra ob- 
tained after subtraction of the pre-burst persistent emis- 
sion from the total source spectrum) by a black body, the 
spectral fits are almost all satisfactory (x 2 e d values of ~1). 
The fit results are shown for the 3 short bursts in Fig. [ll] 
(bl, b3, b5) and for the 7 long bursts in Figs. |l2| (b2, 
b4, b7-bl0) and [l3] (left panel; b6). For the long bursts 
we used a logarithmic time scale, to emphasize the start 
of the bursts, where changes in the parameters are most 
rapid. 

The fits are, while much improved overall, still not op- 
timal during the peaks of the long bursts b4, b6 and blO 
( x 2 ed =2-3 with 44 [b4] or 37 [b6,bl0l dof for the spec- 
tral fits to the 16s spectra). In Fig. [lj we show the aver- 
age net burst spectrum during the flat top part of burst 
b6 excluding the radius expansion and initial contraction 



phase (Sect. 3.2.4), i.e. 27-187s after the start of the burst. 
Tbb and Rbb.io do not change much during this inter- 
val. Clearly, deviations occur below ~10keV; the intensity 
drops below ~5.5keV, while there is an excess between 
~5.5-8keV. These deviations are much larger than the 
calibration uncertainties. 

In the short bursts Tbb decreases during the decay, 
indicating cooling of the neutron star, as already noted 
from the hardness curves. There are some slight variations 
in Rbb.io and Tbb during the first few seconds. However, 
considering the behaviour of the net burst black-body flux, 
Fbb,not, they do not show the correlations that would be 
expected for radius expansion events, i.e. an increase in 
Rbb,io with a simultaneous drop in Tbb and a (nearly) 
constant F b b,not- 

The long bursts all show more or less the same be- 
haviour (note again that for burst b2 we have no short 
time scale spectral information, see Fig. [l2]). After a fast 
increase (^0.5 s) to maximum, Fbb,not remains constant 
(1.4-1.7xl0 -8 ergs -1 cm -2 for the various bursts) for 50- 
200 s, and then decays exponentially. The duration of this 
flat-topped phase is similar to that seen in the light curve 
at high energies (see middle panels of Figs. |] and ||). 
During the first part Tbb increases, then it levels off, to 
decrease again when Fbb,nct decreases. 

Rbb,io is —5 km during the decay of the burst, until 
Fbb,net drops below typically 0.2 xlO -8 ergs -1 cm -2 . Then 
Rbb,io decreases again. Note that this latter behaviour re- 
sembles the systematic effects described by van Paradijs 
& Lewin (1986) near the end of the burst, when the per- 
sistent emission becomes important with respect to the 
net burst emission (Sect. 3.2.2.1). However, we have ar- 
gued above that the burst black-body emission is decou- 
pled from the persistent black-body emission. It is unlikely 
that this behaviour is due to the fact that we assumed a 
constant persistent component during the bursts, while in 
fact the persistent component may vary on the same time 
scale of the burst or faster (see Sect. 3.1). If that was the 
case, one would expect that in some bursts Rbb would de- 
crease, whereas in some others it would increase. On the 
contrary, we find that Rbb decreases in all bursts. We will 
discuss this issue in Sect. 5.2.2. 

3.2.4. Radius expansion bursts 

The most dramatic changes in the spectral parameters can 
be seen during the first few seconds of the long bursts: 
while there is a short dip in Tbb near the start of the 
burst, Rbb, 10 shows values which are up to a factor of ~6 
larger than those found later on (although this is not true 
for burst blO, see Fig. [l2]). This is typical of radius expan- 
sion/contraction episodes, as described in Sect. 1. Burst b6 
shows this behaviour clearest, and may have had even two 
of such episodes within one second (Fig. left panel). 
The expansion to maximum radius takes place in less than 
~0.5s, while the initial contraction phases are longer: it 
takes between 5 and 20 s from maximum expansion back 
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Fig. 12. Spectral fit results for the net burst emission of bursts b2, b4 and b7-bl0 plotted on a logarithmic time 
scale; from top to bottom: bolometric black-body flux, Fbb,nct, in 10~ 8 ergs -1 cm~ 2 , black-body temperature, kTbb, 
apparent black-body radius, Rbb,io, at lOkpc, and goodness of fit expressed in reduced x 2 . The filled circles and 
open squares represent the fit results of the 0.25 s and 16 s spectra, respectively. The data have been logarithmically 
rebinned for clarity. For the 0.25 s spectral fits the number of dof is 18 for burst b4 and 16 for bursts b7-bl0. For the 
16 s spectral fits the number of dof is 44 for bursts b2 and b4, and 37 for bursts b7-bl0. 



to close to the neutron star surface. The final contraction 
phase lasts between ~20 and ~180s We conclude that 
bursts b4, and b6-b9 are radius expansion bursts. 

In Fig. p~o| b we show the (average) net-burst peak 
fluxes, Fbb,net_peak, determined in the same way as in 
the previous subsection for Fbb.pcak, but now for the net- 
burst emission. Most of the net burst peak fluxes are in 
the same range. For the radius expansion bursts this is 
expected, if the Eddington limit is reached. Note that 
Fbb, P cak— Fbb,not_ P cak for bursts b6 and b7, as expected if 
there is no black-body contribution to the persistent emis- 
sion. We conclude that Fbb,not_ pC ak in the radius expansion 
bursts is a measure of the Eddington luminosity observed 
at Earth. We find that Fbb,nct_poak for the non-radius ex- 
pansion bursts (except bursts b3) is more or less similar 
to the maximum flux for the radius expansion bursts. 

Radius expansion events are most conveniently dis- 
played in a flux-temperature diagram (where temper- 
ature increases from right to left, analogous to HR- 
diagrams). Our best example, burst b6, is displayed this 



way in the right-hand panel of Fig. |l3j. In this dia- 
gram the radius expansion/contraction phase and the 
subsequent cooling of the neutron star are clearly dis- 
tinguished by two separate tracks (see e.g. Lewin et al. 
1993). The data points distributed along the horizontal 
line in the upper part of the diagram (i.e. nearly con- 
stant Fbb,nct of 1.2-1.3xl0~ 8 ergs _1 cm~ 2 ) represent the 
expansion/contraction phase. The data points distributed 
along the diagonal line from the upper left to the lower 
right part of the diagram are from the cooling phase of the 
burst. We note that the point in between the two tracks 
is from the very first rise phase. Comparing both pan- 
els of Fig. |l^ it is apparent that the burst spends a long 
time near the vertex of both tracks, about 150 s, where 
Fbb.nct^l^lxlO-Sergs- 1 cm" 2 and kT bb ^2.65 keV. The 
first part (Fbb,nct^0.4x 10~ 8 ergs -1 cm~ 2 ) of the cool- 
ing track is well described (x 2 e d=0-6 for 9 dof) by a 
straight line with a slope of 4.14±0.15 (dotted line in 
Fig. |IJ), i.e. Fbb,nct is consistent with being proportional 
to T^ b , as expected if the neutron star photosphere radi- 
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Fig. 13. Le/£ panel: same as Fig. |l2|, but for burst b6, the longest burst in our sample. Right panel: Bolometric black- 
body flux (Fbb.not in 10 -8 ergs -1 cm -2 ) versus black-body temperature (kTbb) for burst b6. The dotted line represents 
the fit to the cooling track of the burst for data points with Fbb,not > 4 x 10~ 9 ergs -1 cm -2 , see text. Note that kTbb 
increases from right to left. For the 0.25 s and 16 s spectral fits the number of dof is 16 and 37, respectively. 



ates as a black-body with constant effective area. Below 
Fbb.net— 0.4x 10 -8 ergs -1 cm -2 the data start to deviate 
from this relation. The flux-temperature diagrams for 
the other long bursts are more or less consistent with 
the behaviour of burst b6, although the exact locations 
of the vertices between the radius expansion/contraction 
phase and the cooling phase differ slightly (kTbb— 2.3- 
2.6keV, F bb ,nct^l-2-1.35xl0 -8 ergs -1 cm -2 ). The early 
parts of the cooling phases for the other long bursts 
are in most cases also consistent with black-body cool- 
ing, although the flux level at which they start to de- 
viate from Fbb,netCxT^ b varies from burst to burst (e.g. 
~1 x 10 -8 ergs -1 cm -2 for burst b4). The slope becomes 
somewhat steeper below these fluxes for all bursts. 

3.2.5. The radius expansion/contraction track 

As we have shown above, the observed net burst flux 
is consistent with being constant during the radius ex- 
pansion/contraction phases. For instance, for burst b6 a 
constant fit yields Fbb.nct— 1.24xl0 -8 ergs -1 cm -2 with a 



reduced % 2 of 1.4 for 29 dof, excluding the highest flux 
point. However, for a constant composition and constant 
(an)isotropy of the radiation, one would expect to see a 
slight increase in the observed bolometric burst flux with 
increasing photosphere radius, since the gravitational red- 
shift decreases (see e.g. Lewin et al. 1993). In Fig. [l5[ 
we zoom in on the observed radius expansion/contraction 
phase of burst b6 (note that for Fbb,net and Tbb we now 
use a linear scale). Overdrawn is an example of the ex- 
pected relation between bolometric burst flux and the pho- 
tospheric radius (continuous line), for 'standard' values of 
the mass of the neutron, M ns =1.4MQ, the hydrogen frac- 
tion (by mass), A=0.73, the spectral hardening factor, 
Tbb/T e ff,oo=f -7, and high-temperature electron scatter- 
ing opacity, see Appendix B for details. The curve was 
normalized so as to give the observed values of Fbb,net 
and Tbb at touch down (the leftmost point of the radius 
expansion/contraction track). In practice this is analogous 
to solving for the distance and the photospheric radius R 
at touch down (see Appendix B). We can see that the 
curve does not follow most of the data points. Taking the 
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low-temperature electron scattering opacity (but keeping 
the other parameters at their 'standard' values, see above) 
makes the disagreement even larger (dash-dot-dot-dotted 
line in Fig. |l5|). There are three ways to flatten the ra- 
dius expansion/contraction track, or equivalently, to re- 
duce the gravitational redshift corrections. The most ob- 
vious one is to lower the mass of the neutron star; the ex- 
pected track (dash-dotted line) for M ns =0.5MQ is given. 
However, such a low-mass neutron star is in not in line 
with the observed and expected masses for neutron stars 
(see e.g. Thorsett & Chakrabarty 1999). A second option 
is to lower the hydrogen content of the burning mate- 
rial (shown by the dashed line for X=0). However, the 
long contraction phase is typical for unstable mixed H/Hc- 
burning (see Sect. 5.4) and not expected for unstable pure 
He-burning. Moreover, the distance derived at touch down 
is then rather high (d=15.2 kpc for X=0, see Appendix B). 
Finally, increasing the hardening factor to Tbb/T e ff,oo=2 
also flattens the expected track (dotted line in Fig. |l5|) . 
This may be the most realistic option, since spectral hard- 
ening values of ~2 are inferred for burst luminosities near 
the Eddington limit (Babul & Paczyhski 1987; but see 
Titarchuk 1994, where Tbb/T ff,oo<l at very large pho- 
tospheric radii). We note that the effects of an expanded 
boundary layer, as discussed in Sect. 5.2.1, may also be of 
importance here. However, at present it is unclear how to 
take this into account. 

3.2.6. Burst parameters 

From X-ray spectral fits we can determine, for each burst, 
the maximum net bolometric black-body flux, Fbb,maxj 
and the total burst fluence (i.e. the integrated net burst 
flux), Eb- These can be used to derive the burst parame- 
ters 7=F pers /Fbb,max) where F pcrs is the persistent flux 
flux between 2 and 20keV (Table 3), and the average 
burst duration, r=Eb/Fbb,max- Note that Fbb.max may be 
slightly underestimated due to the finite width of the time 
bins. The magnitude of this effect depends on how fast the 
flux varies with time. Eb has been determined by adding 
up the observed net burst fluxes, Fbb,net, per time bin 
from the beginning to the estimated end of the burst. Since 
Fbb,not decays exponentially it only vanishes finishes at in- 
finite times. To compensate for this we fit the decay with 
an exponential and determine the 'rest' fluence by inte- 
grating from the estimated end of the burst up to infinity. 
The 'rest' fluence is large only for burst b2, which was in- 
terrupted during its decay (see Fig. |]) . We also determined 
the burst parameter a, which is the ratio of the average 
persistent flux to the time-averaged flux emitted in the 
bursts, a=F pcrs /(Eb/At), where At is the time since the 
previous burst. For a we can only give lower limits, since 
the source is not observed during South Atlantic Anomaly 
passages and earth occultations. For burst b5 we have, 
however, also assumed that between bursts b4 and b5 no 
other bursts occurred (At=5.77hr). Note that when we 
derive a we assume that between bursts the persistent lu- 



minosity is constant. Since GX17+2 is a highly variable 
source (see e.g. Fig. |l|), this is not strictly valid. The burst 
parameters can be found in Table 4. 

3.3. Burst position in the Z 

The colour-colour diagram (CD) of the data from the three 
different RXTE gain epochs is shown in Fig. [if], together 
with the source positions just before the burst (see also 
Table 2). To determine the position of the source in the 
CD at the time of the burst, we calculated the soft and 
hard colour values from 64 s intervals just before a burst. 
The three short bursts all occurred when the source was 
in the lower part of the NB. Among the long bursts, two 
occurred in the lower part of the HB, four in the NB, and 
one close to the NB/FB vertex. During the observations 
in 1996-2000 GX17+2 spent about 28% of its time in 
the HB, 44% of its time in the NB, and 28% in the FB. 
It thus seems that the NB is overpopulated, whereas the 
FB is underpopulated with bursts (i.e. none occurred). 
However, by assuming that bursts have an equal chance to 
occur at each instant independent of branch, this result is 
statistically not significant (4% probability seeing 7 out of 
10 bursts in the NB, and 7% of seeing none out of 10 bursts 
in the FB). Our conclusion does not change much if we 
include the EXOS AT/ME results reported by Kuulkers et 
al. (1997); for 260ksec of observing time, GX17+2 spent 
11% of the time in the HB, 66% in the NB and 23% in 
the FB, with 3 bursts occurring in the NB and 1 in the 
(lower part of the) FB. This leads to total probabilities 
of 7% of seeing 10 out of 14 bursts in the NB, and 6% of 
seeing only 1 out of 14 bursts in the FB. 

4. Search for burst oscillations 

Following the method described by Leahy et al. (1983; see 
also Vaughan et al. 1994) we defined a 99% confidence level 
above which powers are regarded to be due to a real signal. 
This level depends on the number of trials, which basically 
is the number of independent frequency bins that are ex- 
amined. In our search this number is ~2x 10 7 which leads 
to a trigger level of 42.8 in the Leahy et al. (1983) nor- 
malization. This level was never reached, indicating that 
no coherent oscillations were detected in any of the bursts. 
We inspected the frequencies of the highest observed pow- 
ers. These powers occurred at apparently random frequen- 
cies, consistent with the idea that they are due to random 
fluctuations. 

Since no powers were found above the trigger level, 
only upper limits to the strength of possible burst oscilla- 
tions can be given (see van der Klis 1989; Vaughan et al. 
1994) . We used the (maximum) observed powers to deter- 
mine a 99% confidence upper limit to a true signal power. 
These upper limit powers were then converted into the full 
amplitude of an assumed sinusoidal signal^], which is given 

3 If a sinusoidal signal is given by y = a sin x , the full ampli- 
tude is defined as A = 2a. 
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burst 


TT 1 a 
r bb,max 


E b " 


r (s) 


7 


a 


bl 


1.59±0.06 


5.38±0.07 


3.39±0.14 


1.3±0.2 


£700 


b2 c 


1.33±0.01 


356.4±0.7 


267±2 


1.7±0.4 


£15 


1)3 


0.79±0.05 


5.54±0.09 


7.0±0.5 


2.4±0.4 


£250 


b4 


1.59±0.06 


430.2±1.1 


271±11 


1.2±0.2 


£10 


1)5 


1.46±0.07 


5.26±0.07 


3.60±0.17 


1.6±0.2 


£100 


1)5 










9000±1300 d 


1)6 


1.68±0.09 


664.1±1.9 


396±22 


1.5±0.1 


£1 


b7 


1.36±0.08 


143.3±1.0 


106±6 


1.8±0.1 


£5 


b8 


1.53±0.08 


176.4±1.1 


115±6 


1.6±0.3 


£25 


b9 


1.50±0.08 


137.6±1.0 


92±5 


1.6±0.3 


£25 


blO 


1.65±0.08 


169.0±0.9 


102±5 


1.2±0.2 


£15 



Burst fluence in 10 erg cm - . 
c Using the 16 s spectral fit results. 

d Assuming that between burst b4 and b5 no other bursts occurred. 
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Fig. 14. Top panel: Average observed net burst (i.e. total 
source minus persistent emission) spectrum during the flat 
top of burst b6, i.e. from 27 to 187 s after the start of the 
burst. A black-body fit it shown (kT bb =2.658±0.015keV, 
R bb 10 =4.65±0.06km, x? cd =2.4 for 37 dof). Bottom panel: 
Residuals after subtracting the best fit black-body model 
from the observed spectrum. 
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Fig. 15. Same as the right panel Fig. [j^, but only the 
radius expansion/contraction phase is shown (note that 
the highest flux point lies outside the plot boundaries). 
The five lines correspond to expected tracks in the same 
diagram for different sets of parameters, see text. 
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by: 



A = I.6I1 



where P u \ is the (upper limit) signal power (Vaughan et al. 
1994; Groth 1975), A" the number of photons used in the 
power spectrum, v the frequency at which the power was 
found, and i^sr yq the Nyquist frequency. Table [| gives the 
upper limits on A for all bursts — each burst was divided 
into four parts: 5 s before the rise, the rise, the peak, and 



decay (= end of peak + 2xt oxp ). Note that A is a factor 
\/2 larger than the rms amplitude that is often quoted. 

The values of the upper limits are not very constrain- 
ing. Only those obtained from the 2 s power spectra in 
the 2-60 keV band are well below the strongest values ob- 
served in some of the atoll sources (see e.g. Strohmayer 
1998, 2001). 

5. Discussion 
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Table 5. Upper limits on burst oscillations 



burst 


p 


upper 


limits on fractional amplitude b 


burst 




upper limits on 


fractional amplitude 6 






0.25 s 


2s 


0.25 s 


2s 






0.25 s 


2s 


0.25 s 


2s 






2-60 keV 


2-60 keV 


8-20 keV 


8-20 keV 






2-60 keV 


2-60 keV 


8-20 keV 


8-20 keV 


bl c 


b 


0.28 


0.12 


>1.0 


0.67 


b6 


b 


0.41 


0.14 


0.70 


0.28 




r 


0.23 


0.07 


>1.0 


0.34 




r 


0.28 


0.09 


0.49 


0.24 




P 


0.25 


0.09 


>1.0 


0.45 




P 


0.40 


0.15 


0.63 


0.24 




d 


0.27 


0.11 


>1.0 


0.60 




(1 


0.45 


0.16 


0.89 


0.30 


b2 


b 


0.34 


0.10 


0.66 


0.26 


b7 


b 


0.47 


0.16 


0.78 


0.26 




r 


0.19 


0.07 


0.51 


0.20 




r 


0.30 


0.08 


0.58 


0.15 




P 


0.29 


0.10 


0.55 


0.20 




P 


0.36 


0.13 


0.63 


0.23 




(1 


0.33 


0.11 


0.68 


0.26 




d 


0.46 


0.16 


0.80 


0.34 


b3 


b 


0.36 


0.13 


0.77 


0.25 


b8 


b 


0.41 


0.14 


0.74 


0.26 




r 


0.21 


0.08 


0.51 


0.17 




r 


0.33 


0.10 


0.56 


0.20 




P 


0.32 


0.12 


0.64 


0.24 




P 


0.37 


0.14 


0.59 


0.24 




d 


0.32 


0.12 


0.71 


0.25 




(1 


0.43 


0.16 


0.72 


0.27 


b4 


b 


0.37 


0.14 


0.76 


0.27 


b9 


b 


0.41 


0.13 


0.77 


0.30 




r 


0.30 


0.09 


0.62 


0.23 




r 


0.25 


0.09 


0.56 


0.22 




P 


0.29 


0.11 


0.55 


0.20 




P 


0.38 


0.14 


0.63 


0.22 




d 


0.38 


0.13 


0.75 


0.27 




d 


0.45 


0.17 


0.78 


0.28 


b5 


b 


0.33 


0.11 


0.63 


0.22 


blO 


b 


0.45 


0.19 


0.67 


0.34 




r 


0.24 


0.09 


0.42 


0.17 




r 


0.34 


0.09 


0.48 


0.23 




P 


0.27 


0.10 


0.39 


0.17 




P 


0.38 


0.14 


0.51 


0.26 




d 


0.28 


0.10 


0.57 


0.22 




d 


0.48 


0.18 


0.68 


0.34 



a Phase of burst profile: b=before burst, r=rise, p=peak, d=decay. 
6 99% confidence upper limits in the 50-2000 Hz frequency range. 

c 99% confidence upper limits in the 2-60 keV energy band are between 50-250 Hz. The high-energy band used is 13.5-20 keV. 



5.1. X-ray burst properties 

Of the ten X-ray bursts from GX17+2 seen so far with 
the RXTE/PCA, three are rather short (~10s), whereas 
the other seven have durations of ~6-25min. They all 
show evidence for cooling of the neutron star photosphere 
during the decay, confirming previous results. The burst 
durations found are in the range seen previously (Kahn & 
Grindlay 1984; Tawara et al. 1984c; Sztajno ct al. 1986; 
Kuulkers ct al. 1997). Using the time-resolved X-ray spec- 
tral fits (see Sect. 5.2) we determined the burst fluence, 
r (a characterisation of the burst decay time), 7 (the net 
brightness of the burst with respect to the persistent level) 
and a (the ratio between the total net burst flux and to- 
tal persistent flux in between bursts). Since not all these 
parameters were given for GX17+2 bursts observed pre- 
viously, we (re-)determined these values from the infor- 
mation given in the various papers. This is explained in 
Appendix A and the results are given in Table A.l. Again, 
our values for the burst parameters for both the short and 
long bursts are comparable to those for the bursts ob- 
served previously. Since observations with RXTE are in- 
terrupted by either SAA passages or Earth occultations, 
and the fact that the length of each 'continuous' observa- 
tion was limited (~hours) we could not determine the time 
interval between bursts. Our lower limits on a are, there- 
fore, rather low. However, from the bursts observed previ- 
ously, it is found that a always exceeds ^1000. This means 
that most of the accreted material is being burned con- 
tinuously. One of our short bursts occurred about 5.8 hr 



after a long burst. So far this is the shortest recurrence 
time encountered for a burst in GX 17+2. The fluences of 
the bursts from GX17+2 are either ~5.5xl0~ 8 crgcm~ 2 
(short bursts) or 140-670xl0 -8 ergcm -2 (long bursts). At 
a distance of lOkpc (but see Sect. 5.3) these correspond 
to total energies of 6.5 x 10 38 erg and 2-8 x 10 40 erg, respec- 
tively, in the range seen for type I bursts (see e.g. Lewin 
et al. 1993). 

One of the short bursts was rather faint, with a peak 
intensity which was a factor of ~1.25 lower than for the 
other bursts, and displayed two peaks in the light and 
bolomctric flux curves. This burst resembles the (rela- 
tively) weak double-peaked bursts seen from 4U 1636— 53 
(Sztajno et al. 1985; Fujimoto et al. 1988). Such bursts 
are clearly not radius expansion/contraction events, but 
are instead thought to be due to the mixing of fresh ma- 
terial into the unstably burning layer; such mixing may be 
caused by shear instabilities in the outer envelope of the 
neutron star (Fujimoto et al. 1988). 

We show for the first time that almost all of the long 
bursts of GX 17+2 have episodes of radius expansion and 
contraction of the neutron star photosphere. In all cases 
the expansion and initial contraction occur during the first 
few seconds of the bursts, but the final contraction phase 
duration varies between ~20s and ~180s. During the ex- 
pansion/contraction phase the net burst flux remains con- 
stant, presumably at the Eddington value. Long contrac- 
tion phases might be caused by long phases of unstable 
hydrogen burning, whose time scale is limited by /3-decays, 
up to ~100s after the burst has started (see e.g. Bildsten 
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Fig. 16. Colour-colour diagrams of the RXTE/PCA observations of GX17+2 (after Homan et al. 2001). The left, 
middle and right panels refer to the observations done during RXTE gain epochs 1, 3 and 4, respectively. The soft 
colour is defined as the count rate ratio in the 4.8-7.2 keV to the 3.0-4.8 keV energy bands for the epoch 1 data, in 
the 4.8-7.3 keV to the 3.0-4.8 keV energy bands for the epoch 3 data, and in the 4.6-7.1 keV to the 3.0-4.6 keV energy 
bands for the epoch 4 data. The hard colour is defined as the count rate ratio in the 10.6-19.7keV to the 7.2-10.6keV 
energy bands for the epoch 1 data, in the 10.5-19.6 keV to the 7.3-10.5 keV energy bands for the epoch 3 data, and 
in the 10.5-19.5 keV to the 7.1-10.5 keV energy bands for the epoch 4 data. Each point represents an average of 16 s. 
We have indicated the position of the source during 64s intervals just before a burst (indicated by bl, b2, blO, in 
chronological order), and the three limbs of the Z (HB, NB, FB). Since the HB is not clearly distinguishable from the 
NB in this figure, we have indicated the approximate HB/NB vertex by a line. 



1998, 2000, and references therein). The maximum ap- 
parent radius at lOkpc, Rbb,i0j reached during the ex- 
pansion episodes is ~30km. Corrected for gravitational 
rcdshift effects and spectral hardening (see Appendix B), 
the maximum photospheric radius R (as measured by a 
local observer) is ~90-125km. We do not see evidence 
for a change in the cut-off power-law component above 
^10keV when the photosphere expands up to these radii. 
This suggests that either the environment in between the 
neutron star surface and expanded photosphere at ~90- 
125 km is involved in producing the high-energy radiation 
but was not affected, or the inner disk and corona which 
presumably produce the high-energy radiation do not ex- 
tend down to the neutron star. A possibility is that the 
magnetospheric radius extends out to at least these radii, 
disrupting the disk and corona near that point. Note that 
the radius range of ~90-125 km overlaps with the range of 
magnetospheric radii expected for Z sources (~ 10-100 km, 
see Ghosh & Lamb 1992). However, at the high mass ac- 
cretion rates inferred for GX17+2 the disk is expected 
to extend down to close to the neutron star (see Popham 
& Sunyaev 2001). For a more detailed discussion of the 
different emission regions we refer to the next subsection. 

Long (>minutes) X-ray bursts have been seen in 
many other sources (see, in order of decreasing dura- 
tion, Hoffman et al. 1978b, Lewin et al. 1984, ~25min: 
4U 1708-23; Swank et al. 1977, ^lOmin: 4U 1724-307; 
Tawara et al. 1984a,b, ^4.5 min: 3A 1715-321: Kaptein et 
al. 2000, ^3.3 min: 1RXS J171824.2-402934; van Paradijs 
et al. 1990, ^2.7 min: 4U 2129+11). The shapes of the light 



curves of these bursts at different energies are very similar 
to the long bursts of GX 17+2, i.e. fast (several seconds) 
rise and peaked emission at low energies and relatively 
slower rise at high energies. They also have rather long 
episodes of contraction, i.e. long episodes of Eddington 
fluxes, up to ~300s (see e.g. Lewin et al. 1984). Moreover, 
the curves for Tbb and Rbb as a function of time show sim- 
ilar shapes. 

The difference, however, between the long bursts of 
GX17+2 and those in the other sources is that 'precur- 
sors' are present in the light curves before almost all of 
the long bursts from other sources. Note that these are 
not real precursors, but indicate very large radius expan- 
sion (Rbb > 100 km). No emission is then seen because of 
the very low effective temperatures reached; the peak of 
the emission is shifted to UV wavelengths. Thus the X- 
ray burst light curve is interrupted by a gap of no X-ray 
emission leading to an apparent precursor (see Lewin et 
al. 1984). The rise to maximum radius in those sources 
is on the order of several seconds. In GX 17+2 the maxi- 
mum value of Rbb is a few times the neutron star radius 
and the expansion lasts typically only a second. The most 
important difference, however, is that the persistent emis- 
sion just before the long bursts in other sources is ~1% of 
the Eddington value, i.e. a factor of ^100 lower than in 
GX17+2 (see also Sect. 5.4). 

We note that even longer X-ray bursts do exist, 
which last for several hours (Cornelisse et al. 2000, 2001; 
Wijnands 2001; Strohmayer & Brown 2001; Kuulkers 
2001; Kuulkers et al. 2001). Their recurrence time is much 
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longer than normal type I bursts (<;7.5 days, Cornclissc et 
al. 2001; ^4.8yrs, Wijnands 2001), and they seem to occur 
at persistent luminosites which are a factor of ~0. 1-0.2 of 
the Eddington limit. 

5.2. Black- body spectral fits 

5.2.1. Decoupled persistent and burst emission 

It has previously been pointed out (van Paradijs & Lewin 
1986; see also Sztajno et al. 1986) that, if the persistent 
emission before the burst contains a spectral component 
which originates from the same regions as the burst emis- 
sion, this should be taken into account in the spectral 
analysis of X-ray bursts. This is especially the case when 
the neutron star is accreting matter at high rates such 
as thought to be the case in Z sources, making the neu- 
tron star hot and radiate soft (black-body like) emission. 
'Standard' X-ray burst spectral analyses (modeling the 
net-burst emission by black-body radiation) then results 
in systematic errors in the spectral parameters, especially 
during the end of the burst when the net-burst luminos- 
ity is small. It was therefore suggested to model the total 
burst emission using a two-component approach, i.e. with 
a soft component (usually a black body) and a hard com- 
ponent (fixed at the parameters found for the persistent 
emission). The former component is thought to contain 
all emission from the neutron star photosphere, while the 
latter component is assumed to be from the accretion disk 
and not affected by the prompt burst. However, we find 
that using this two-component spectral fitting approach 
does not give satisfactory results. Instead, the 'standard' 
burst analysis, which assumes the persistent emission con- 
tinues unchanged, yield better fits. 

This suggests that the persistent emission is unaffected 
during the burst and that the persistent black-body com- 
ponent (which is not always present) does not arise from 
the same region as the burst emission, contrary to what 
was previously assumed. We find that during the initial 
exponential decay phase of the bursts the neutron star 
photosphere shows pure black-body cooling. During that 
stage the apparent black-body radius stays constant at 5- 
6 km (at lOkpc). Comparing these radii with the apparent 
radii of the black-body component of the persistent emis- 
sion (if present) indicates that the persistent black-body 
component originates from a larger area than that of the 
burst emission (see Fig. ||). The inner part of the accretion 
disk or an expanded boundary layer are the most obvious 
sites for this persistent soft component. A recent study 
of X-ray spectra from LMXBs with luminosities spanning 
several orders of magnitude has shown that the persis- 
tent black-body emission likely does not arise in (part of) 
the accretion disk (Church & Balucihska-Church 2001). 
Instead at very high accretion rates (near Eddington val- 
ues) the accretion disk boundary layer around the neutron 
star expands radially, with radial extents larger than one 
stellar radius (Popham & Sunyaev 2001). Moreover, at 
these high rates the expected spectrum of the boundary 



layer strongly resembles a black-body. The fact that we see 
black-body emission in the persistent spectra with appar- 
ent radii which are larger (factor ^2-3) than the neutron 
star (see also Church & Balucihska-Church 2001) is in 
line with this. Still, one would expect that such a bound- 
ary layer would be affected by the bursts (see Popham 
& Sunyaev 2001), of which we do not see clear evidence 
from our spectral fits. We note, however, that Homan et 
al. (2001), who studied the properties of the NB and FB 
quasi-periodic oscillations (NB/FB QPO or NBO/FBO) 
during bursts b4 and blO, found that during the bursts the 
absolute amplitude of these QPOs decreased significantly. 
This suggests that the inner accretion flow is affected by 
the increase in the radiation from the neutron star, con- 
trary to what our spectral results indicate. Apparently the 
QPO mechanism of the NBO/FBO is much more sensitive 
to the radiation field than the bulk flow of matter itself. 
This is in accordance with models for NBO /FBO (Fortner 
et al. 1989), which require a delicate balance between the 
radiation field and the accretion flow, that might easily be 
disrupted during X-ray bursts. 

5.2.2. Deviations from black-body emission 

During the radius expansion/contraction phase black- 
body emission does not provide a good description of 
our observed net-burst spectra, leaving a slight excess 
in emission between ~5.5-10keV and lack of emission 
between ~3-5.5keV (Fig. 11). Such systematic devia- 
tions have been observed before, especially during radius 
expansion/contraction phases. For example, spectra ob- 
tained with the Large Area Counter onboard Ginga dur- 
ing the long burst of 4U 2129+11 (van Paradijs et al. 1990) 
show remarkably similar deviations at the same energies 
as we see in GX17+2. The persistent mass accretion in 
4U 2129+11 is inferred to be a factor of ~100 lower than 
for GX 17+2. This suggests that the radiation properties 
of the photospheres are the same when they are radiating 
near the Eddington luminosity, which is independent of 
the mass accretion rate. As noted by van Paradijs et al. 
(1990), the 'bumpiness' of the residuals point to the pres- 
ence of relatively narrow-band spectral features, rather 
than the residuals being due to broad-band spectral hard- 
ening at Eddington luminosities (see also Strohmayer & 
Brown 2001). 

From our 'standard' spectral analysis, we find that 
the apparent black-body radii decrease near the end of 
the decay of the bursts. This was also inferred from the 
EXOSAT/ME observations of GX17+2 (Sztajno et al. 
1986), and has been seen in other sources as well (e.g. 
Chevalier & Ilovaisky 1990). If this were due to a real de- 
crease in the burst emitting area one might expect to see 
burst oscillations due to the spinning neutron star. We do 
not see evidence for this. If the expanded boundary layer 
reacts to its changing environment on the same time scale 
as the burst itself then variations in this layer's structure 
could also play a role here. 
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This observed behaviour resembles the systematic ef- 
fects described by van Paradijs & Lewin (1986), which 
occurs near the end of the bursts in the presence of a 
persistent hot neutron star component. However, we have 
argued that the persistent black body component does 
not contribute to the burst emission, so this is not a vi- 
able explanation. Sztajno et al. (1986) argued that this 
effect also can not be due to spectral hardening, since that 
would lead to an apparent radius increase with decreasing 
Tbb- Their argument was based on the assumption that 
the spectral hardening decreases with decreasing tempera- 
tures (and therefore decreasing burst luminosity; see also 
Sztajno et al. 1985). However, later model calculations 
have shown that in fact spectral hardening increases again 
whenever the burst luminosity drops below a certain value 
(L b b^0.2L E dd, e.g. London et al. 1986, Ebisuzaki 1987). 
For an assumed constant size of the emission region dur- 
ing the cooling stage (presumably the whole neutron star 
surface), the apparent radii will then become smaller as 
Tbb decreases, consistent with what is observed. 

5.3. The distance to GX 17+2 and the persistent 
emission 

Assuming that during the expansion and contraction 
phase the net-burst luminosity equals the Eddington lu- 
minosity we can estimate the distance to GX 17+2. As 
noted before (Lewin et al. 1993; see also recent discus- 
sion by Kuulkers & van der Klis 2000), this is not with- 
out problems. Assuming (see Appendix B) a neutron 
star with a canonical mass of 1.4 Mq, using the elec- 
tron scattering opacity, k, for cosmic composition in the 
low-temperature limit, we obtain LEdd = 4:TtcGM/k = 
2 x 10 38 ergs _1 . For the observed values of F bo i— 1-2- 

I. 3x 10~ 8 ergs -1 cm~ 2 this leads to distances of d~11.4- 

II. 9 kpc (assuming isotropic radiation and no spectral 
hardening, i.e. T e ff j00 =Tbb)- Close to the neutron star, 
however, one still has to take into account gravitational 
rcdshift effects. We therefore estimated the distance by 
solving the relevant equations, for the explanation of 
which we refer the reader to Appendix B. Correcting 
for gravitational redshift effects, taking the same 'stan- 
dard' parameters as above, using values of n appropriate 
for very high temperatures, and correcting for spectral 
hardening (Tbb/T e ff )0 o=l-7) we found distances between 
10.8-12.0kpc, with an average of d~11.3kpc (again see 
Appendix B for details). If the burst emission is highly 
anisotropic (£—2) we derive d~8kpc, see below. We note 
that these estimates have systematic uncertainties of the 
order of 30% or so (see Lewin et al. 1993; see also Kuulkers 
& van der Klis 2000, and Appendix B). At touch down the 
radius of the photosphere, R, presumably equals the radius 
of the neutron star surface. We find values of R~ 12 20 km, 
where the range arises from the different assumptions in 
deriving these radii (Appendix B). 

One can also use the radius expansion events to esti- 
mate the persistent luminosity in Eddington units, with- 



out the need of knowing the distance to the source, mass 
of the neutron star, etc. However, for the long bursts we 
find that the total emission before/after the burst is a up 
to factor of ~2 larger than the net-burst emission during 
expansion/contraction phase. If indeed the maximum net- 
burst flux equals the Eddington flux, this implies that the 
persistent flux is up to 2 times the Eddington flux. From 
the models characterising the behaviour of Z sources it 
is inferred that on the HB, NB and FB, the mass accre- 
tion rates (and therefore luminosity) are just below, at or 
just above Eddington values, respectively (e.g. Hasinger 
1987; Lamb 1989; Hasinger et al. 1990). The only way 
out is that the net burst emission is (highly) anisotropic 
compared to the persistent emission. Since the true lumi- 
nosity values should be more or less comparable, one thus 
infers an anisotropy factor of ^2 for the net burst emis- 
sion. This cannot be ascribed to localized emission on the 
neutron star (see also Section 5.4), but might be due to 
part of the burst emission being hidden behind a puffed- 
up inner accretion disk. In the case of highly anisotropic 
burst emission (£—2) the distance estimate to GX17+2 
is reduced to ~8kpc (see Appendix B). This estimate is 
comparable to that derived for other similar X-ray sources 
(e.g. Christian & Swank 1997; Schulz 1999). We note that 
anisotropic burst emission may also explain the apparent 
inconsistency in the distance to Cyg X-2 from a radius ex- 
pansion burst, i.e. ~11.6kpc (Smale 1998; for a 1.9 Mq 
neutron star, peak photospheric radius of 26 km, cosmic 
abundances, low temperature opacity, no spectral hard- 
ening; see also Appendix B), and that derived from opti- 
cal measurements, i.e. ~7.2kpc (Orosz & Kuulkers 1999). 
However, for Cyg X-2 the persistent flux before the burst 
was of the same order as the net peak burst flux. 

5.4. X-ray bursts and (extreme) mass accretion rates 

We found no significant correlation between the burst pa- 
rameters and the position of the source in the Z-track at 
the time of the bursts, i.e. both long and short bursts 
may occur at similar mass accretion rates. We find that 
no X-ray bursts occurred when the source was on the 
FB, despite considerable coverage compared to the other 
branches of the Z. However, this result is not sufficiently 
significant to exclude the possibility that bursts occur on 
the FB as frequently as elsewhere. We do note that the 
absence of bursts on the FB is consistent with theoreti- 
cal expectations (see Sect. 1; e.g. Bildsten 1998, 2000), as 
the mass accretion rate on the FB is thought to be super- 
Eddington (e.g. Hasinger 1987; Lamb 1989; Hasinger et 
al. 1990). 

At ncar-Eddington accretion rates, long bursts are ex- 
pected from X-ray burst theory (see Sect. 1; see also van 
Paradijs et al. 1988). These are thought to be due to mixed 
H/Hc burning, triggered by thermally unstable He igni- 
tion. The fast rise and the short radius expansion part of 
the long bursts resemble pure helium flashes which last on 
the order of 5-10 s, although in the case of GX 17+2 the re- 
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currence times of these bursts are on the order of a day in- 
stead of the expected hours for typical pure helium flashes. 
This flash triggers the long phase of unstable mixed H/He 
burning. Our analysis provides the first example of the ex- 
istence of such long bursts at high (near-Eddington) accre- 
tion rates. However, GX17+2 also displays short (~10s) 
bursts, which are not expected at high accretion rates. 
Similarly, CygX-2, also accreting at near-Eddington val- 
ues, displays rather short bursts (~5s). Such short bursts 
(presumably He flashes) are only expected for accretion 
rates which are a factor of 20-100 lower. 

Since it is not the global mass accretion rate that 
matters, but the mass accretion rate per unit area (e.g. 
Marshall 1982; see also Bildsten 2000), it might be that at 
different times the area which accumulates most of the ac- 
creted matter differs, giving rise to short and long bursts. 
However, if the areas are relatively small, pulsations are 
expected due to the neutron star rotation (unless the ac- 
creted matter is distributed symmetrically along the ro- 
tation axis, such as in an equatorial belt). We do not see 
any evidence for this. 

The properties of X-ray bursts depend not only on the 
mass accretion rate (per unit area), but also on the tem- 
perature of the neutron star envelope and composition of 
the accreted material (e.g. Fushiki & Lamb 1987; Taam 
et al. 1996; see also Bildsten 1998, 2000; Lamb 2000). 
At high accretion rates and high envelope temperatures 
the combined He flash and mixed H/He burning occurs; 
however at high accretion rates and low envelope temper- 
atures only He flashes may occur (see Lamb 2000). The 
difference in envelope temperature at different times may 
indeed explain the occurrence of the two types of bursts, 
and it would be worthwhile to explore this further. It is 
not clear, however, how such different envelope temper- 
atures could be reached at very similar accretion rates. 
The composition of the accreted material is not expected 
to change much, since fresh matter arrives at high rates. 

The situation at high mass accretion rates becomes 
even more confusing when one considers the bright 'GX 
atoll' sources (GX3+1, GX13+1, GX9+1, GX9+9) as 
well as the other Z sources (ScoX-1, GX5-1, GX 340+0, 
GX 349+2, CygX-2). The former are thought to accrete 
with rates around 10% of the Eddington mass accretion 
rate (Psaltis & Lamb 1998; see also Kuulkers & van der 
Klis 2000). For those sources, only short (~10-15s dura- 
tion) bursts are seen (infrequently) in GX 3+1 (Makishima 
et al. 1983; Asai et al. 1993; Molkov et al. 1999; Kuulkers & 
van der Klis 2000) and GX 13+1 (Matsuba et al. 1995) and 
no bursts have been reported for GX 9+1 and GX 9+9, de- 
spite ample observing times. 

The X-ray spectral and fast timing properties of the 
other five canonical Z sources are very similar to those 
of GX17+2 (see Hasinger & van der Klis 1989; van der 
Klis 2000), and thus one would infer comparable mass ac- 
cretion rates. Of the five of these sources, however, only 
CygX-2 shows X-ray bursts, which are typical He flashes, 
whereas the others do not burst at all (despite ample ob- 
serving times). Since at very high accretion rates no bursts 



are expected one would then naively infer that the aver- 
age mass accretion rate is different, i.e. higher, in the four 
non-bursting Z sources with respect to the two bursting 
Z sources. But this would imply that the correlated X-ray 
spectral and fast timing properties are not a function of 
the inferred mass accretion rate; van der Klis (2001) has 
recently proposed a way in which the correlated properties 
could vary in response to changes in the mass accretion 
rate without being a function of it. A slightly different 
suggestion was made by Homan et al. (2001) who propose 
that the mass accretion rate is not the parameter that de- 
termines the position along the Z track. As a consequense, 
bursts that occur at a similar position along the Z might 
occur at different mass accretion rates and therefore have 
different proporties. Note also that differences in mass ac- 
cretion rate per unit area may be of importance here, as 
was already suggested by Kuulkers et al. (1997) to explain 
the properties of quasi-periodic oscillations occurring on 
the HB and upper parts of the NB, and the presence of 
bursts in GX 17+2. On the other hand, since the envelope 
temperatures and composition also influence the bursting 
properties, they may be of importance too (e.g. Taam et 
al. 1996). Future modelling may provide more insight in 
the processes involved. 

6. Conclusions 

• We found ten X-ray bursts in all data on GX 17+2 ob- 
tained with the RXTE/PCA to date. Three of them were 
short (~10s), while the others lasted from ~6 to 25min. 
All the bursts showed spectral softening during the de- 
cay, indicative of cooling of the neutron star. No evidence 
for high-frequency (>100Hz) oscillations at any phase of 
the bursts is seen. We find no evidence for correlations 
of the burst properties with respect to the source posi- 
tion in the colour-colour diagram (presumably a function 
of the mass accretion rate). The long bursts are consis- 
tent with being due to mixed H/He burning, triggered 
by thermally unstable He ignition, as expected at the in- 
ferred near-Eddington mass accretion rates in GX 17+2. 
However, the presence of short bursts in GX 17+2, as well 
as short bursts in CygX-2, another persistently bright 
LMXB, and no bursts in the other four similar LMXBs 
ScoX-1, GX5-1, GX340+0 and GX349+2, is not ac- 
counted for in the current X-ray burst theories. Note that 
this also holds for the bright 'GX atoll' sources GX3+1, 
GX13+1, GX9+1 and GX9+9, which arc thought to 
accrete near one tenth of the Eddington rate. Of these 
four atoll sources, only GX 3+1 and GX 13+1 infrequently 
show short (~10-15s) X-ray bursts, whereas GX9+1 and 
GX9+9 have not been seen to burst at all. 

• We found that two-component spectral fits to the 
total burst emission, as has been suggested previously, do 
not give satisfactory results whenever the persistent emis- 
sion before the burst contains a black-body component. 
On the other hand, the 'standard' spectral fit analysis, in 
which the burst emission is modeled by black-body radi- 
ation after subtraction of the persistent emission before 
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the burst, does provide satisfactory results. This means 
that whenever there is a black-body contribution in the 
persistent emission before the burst, it is probably also 
present during the burst itself (i.e. the burst black-body 
emission appears on top of the persistent black-body emis- 
sion). This implies that the burst emission does not arise 
at the same site as the black-body component often seen in 
the persistent emission. We find evidence that the emission 
region of the persistent black-body component is larger 
than that of the burst emission, which indicates that it 
probably originates in an expanded boundary layer, as re- 
cent modeling indicates. 

• Five of the long bursts showed evidence for an 
expanding photosphere during the first seconds of the 
burst, presumably due to the burst luminosity reach- 
ing the Eddington limit. The contraction phase differs in 
duration between the bursts, the longest being ~3 min. 
The net burst fluxes reached during the radius expan- 
sion/contraction phase are all more or less similar, which 
is indicative of the Eddington limit as observed at earth. 

• The persistent flux just before the burst is inferred 
to be up to a factor 2 higher than the net burst flux 
during the radius expansion/contraction phase. If in both 
cases the emission is isotropic, this implies that the per- 
sistent emission is up to a factor of 2 higher than the 
Eddington luminosity. We suggest that the burst emis- 
sion is (highly) anisotropic. Since both the persistent lu- 
minosity in GX 17+2 and the net burst luminosity during 
the expansion/contraction phase are thought to have near 
Eddington values, we then derive an anisotropy factor of 
—2 for the burst emission. 

• When the burst luminosity is near Eddington val- 
ues, deviations from pure black-body radiation are evi- 
dent below lOkeV. Similar deviations have been seen dur- 
ing (long) X-ray bursts in other low-mass X-ray binaries, 
and can not be explained by spectral hardening. Assuming 
that the black-body approximation is, nevertheless valid, 
that the net burst peak fluxes during the radius expan- 
sion/contraction phase equals the Eddington limit as seen 
on Earth, and using an anisotropy factor of 2, we estimate 
the distance to GX 17+2, taking into account gravitational 
redshift effects and spectral hardening. For 'standard' pa- 
rameters (M ns =1.4MQ, cosmic composition) we derive a 
distance ~8kpc, with a systematic uncertainty of up to 
-30%. 
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Appendix A: X-ray bursts from GX 17+2 observed 
with other instruments 

X-ray bursts from GX17+2 have been previously seen 
with the Monitor Proportional Counter (MPC) onboard 
the Einstein Observatory (Kahn & Grindlay 1984), the 
second Fine Modulation Collimator (FMC-2) onboard 
Hakucho (Tawara et al. 1984c) and the Medium Energy 
(ME) experiment onboard EXOSAT (Sztajno et al. 1986; 
Kuulkers et al. 1997). Not all the bursts parameters as tab- 
ulated in Table 4 for the bursts seen with the RXTE/PCA 
have been given by these authors. In order to compare our 
bursts in more detail with the other observed bursts we 
calculated these parameters from the information given in 
those papers. Since the X-ray energy ranges in which the 
different instruments operate are comparable, and most 
of the (burst) emission is radiated in these energy bands, 
such a comparison is feasible. In the next subsections we 
describe how the parameters were determined for the dif- 
ferent bursts. 

A.l. MPC/ Einstein 

Kahn & Grindlay (1984, see also Kahn et al. 1981) ob- 
served a burst from GX17+2 in 1980 with the MPC on- 
board Einstein, with a decay time of ~7.5s. From their 
background subtracted count rate profiles in the 1.4- 
14.4 keV energy band (their figure 1) we determined the 
values for 7 and r. The value of 7 is consistent with 
what is inferred using the quoted net peak and persis- 
tent count rate of 0.15 and 0.56 Crab (2-6 keV), respec- 
tively, by Kahn et al. (1981). The persistent flux be- 
fore the burst was given as 1.56±0.02xl0 _8 ergcm~ 2 s _1 
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Table A.l. Burst parameters of other bursts in the literature 



blirst Fpcrs Fvjk^ax 


E b c 


r(s) 


7 


a 


Kahn & Grindlay (1984) 










1.56±0.02 see text 


see text 


6.7±1.4 


3.4±0.2 


£1750 


Tawara et al. (1984c) 










A ~2.8 ~1.5 


~346 


367±43 


1.83±0.19 


~1000 


B ~2.8 ~2.0 


~94 


75±8 


1.38±0.13 


~1300 


C ~2.9 ~1.7 


~114 


109±37 


1.77±0.20 


~1000 


D ~2.9 ~2.1 


~406 


309±27 


1.41±0.11 


~1000 


Sztajno et al. (1986) 










1 1.51±0.13 1.1±0.4 


4.7±1.0 


6.27±1.2 


3.6±0.2 


£1450 


II 1.86±0.15 0.99±0.05 


114±3 


132±4 


2.13±0.02 


£650 


Kuulkers et al. (1997) 










III 




~27 


~2.7 


£7000 


IV 




~100 


~3.2 


~2200 


a Unabsorbed persistent flux in 


10" s ergs" 


_i cm -2 , see 


also text. 




6 Peak net black-body flux in 1C 


_s ergs -1 


cm -2 . 







c Burst fluence in 10 erg cm 



(l-20keV). Kahn & Grindlay (1984) report a maximum 
peak flux of 9.1x 1(T 8 erg cm" 2 s _1 (1-20 keV). No er- 
ror is given, but only a ±-sign. It is not clear whether 
this is a typo, or if this means that the value was ap- 
proximate. Their peak flux is much higher than observed 
for the other bursts. An averaged net spectrum over the 
burst profile (first 7.68 s) was modeled as black-body emis- 
sion. Values of 3.8 ^ Tbb ^ 9.3 keV and black-body 
emission area, Abb, of 0.56 ^ Abb ^ 3.2xl0 12 cm 2 (at 
lOkpc) were found. This, in principle, may give us a 
handle on the observed fluence. However, this gives us a 
value of E b of ~161x 10 _8 ergcm -2 (using T b b— 6.55keV 
and Abb— l-88x 10 12 cm 2 ), which is also much larger than 
observed for the other short bursts. We, therefore, do 
not quote the peak flux and burst fluence in Table A.l. 
Using the Einstein database maintained at HEASARCQ 
we found the start of the (simultaneous HRI) observation 
as being March 29, 1980, 00:48:06 (UTC). The burst was 
observed ~57min later (Kahn et al. 1981). Using this and 
the observed rate just before the burst we determined an 
upper limit on a. 

A.2. FMC-2/Hakucho 

Four bursts (two in 1981 and two in 1982; denoted A-D) 
were found, originating from GX 17+2, during pointed ob- 
servations with the FMC-2 onboard Hakucho (Tawara et 
al. 1984c). They had e-folding decay times of ~100s (B, 
C) or -300 s (A, D). For each burst Tawara et al. (1984c) 
tabulated the net peak, net total (=fluence) and persis- 
tent flux in FMC-2 ctss" 1 in the 1-22 keV band. From 
this they calculated 7. Using this information we calcu- 
lated t and verified 7. Tawara et al. (1984c) also provided 
the total integrated energy fluxes (= total fluence) of the 
two bursts in 1981 and the two bursts in 1982, together 
with the integrated persistent fluxes over 40 hr and 50 hr, 
respectively. Note that it is not clear if the fluence and per- 



See http://heasarc.gsfc. nasa. gov / 



sistent flux are both given for the 1-22 keV energy band, 
and whether the persistent fluxes are corrected for inter- 
stellar absorption or not. From this it was estimated that 
a was of the order of 1000. Since the net fluence in counts 
of all bursts individually are given, we can determine the 
individual fluence in energy flux. Also, the (mean) persis- 
tent flux during the two observations can be determined. 
Moreover, assuming that between bursts A and B, which 
occurred 12.5 hr after each other, no other bursts occur, 
we give a more exact value for a in this case (assuming 
a constant persistent flux). For the others we just assume 
the given value of 1000. We also give an estimate of the 
net peak burst flux, from the observed values of F pers and 
7- 



A. 3. ME/EXOSAT 

Two bursts were found in the EXOSAT/ME observations 
taken in 1984 and 1985 with durations of ~10 s and >5 min 
(Sztajno et al. 1986; denoted as burst I and II by Kuulkers 
et al. 1997). The persistent flux (2-20 keV) before a burst 
was determined from the spectral fits done by Sztajno et 
al. (1986), who used a bremsstrahlung plus a black-body 
component, subjected to interstellar absorption. The val- 
ues for the bremsstrahlung parameters and Nh were given 
in their text, whereas the black-body results were indi- 
cated in their figures 4 and 5. The spectral fit results to the 
net burst emission are shown in their figures 2 and 3, from 
which we determined the bolometric maximum net peak 
flux and the fluence. To determine the 'rest' fluence, we 
used the same approach as mentioned in Sect. 3.2.5, and 
fitted an exponential to the burst flux decay. In the case of 
burst I we have only 3 values; we therefore assumed that at 
infinity the burst flux decays to zero, while the exponen- 
tial start time was set to the start of the burst. Since the 
spectral time resolution is not sufficient, we determined 
the burst parameters 7 and r from the observed raw (i.e. 
not dead-time corrected) count rate profiles (1-20 keV) 
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given in their figure 1. We corrected the light curves for the 
'variable' dead time (van der Klis 1989), as appropriate for 
the EXOSAT High Energy Resolution (HER) modes (see 
e.g. Kuulkers 1995). Sztajno ct al. (1986) noted that in 
both bursts the maximum flux is ~40% above the persis- 
tent level, indicating 7~2.5, in rough agreement with our 
estimates. Using the start time of the two observations 
during which the bursts occurred as given by Kuulkers 
et al. (1997), and assuming that the persistent flux was 
constant, we were able to determine lower limits on a. 
Note that Sztajno et al. (1986) gave rough upper limits of 
~10000 and ~500, for burst I and II, respectively, which 
are more or less consistent with our determined values. 

For the two bursts (denoted as burst III and IV) ob- 
served in 1986 discussed by Kuulkers et al. (1997) we only 
estimated 7 and r from the observed (dead-time and back- 
ground corrected) count rate profiles in the 1-20 keV band 
as displayed in their figure 7. r is determined using the in- 
tegrated net burst count rate profile and peak count rate. 
No spectral information is available for the persistent or 
burst emission in order to determine Fbb,max, Eb and a, 
although they note that a is probably larger than ^1000 
for the second burst, if no bursts occurred in between the 
two bursts. However, if we naively assume that the per- 
sistent count rate is constant between bursts, and that 
between the start of the 1986 EXOSAT observation and 
the first burst, and between the two bursts, no other burst 
occurred, we can still get a crude estimate of (the lower 
limit on) the value of a. 



Appendix B: Distance estimate and gravitational 
redshift effect 

Most of the distances derived from radius expansion bursts 
in the literature are solely based on the observed fluxes at 
maximum expansion, Fbb,max, and assuming this equals 
the (non-relativistic) Eddington luminosity, L E dd,non_rcb 
from a source at distance d. Here one uses LEdd,non_rci = 
4TTcGM ns /n, where c is the speed of light, G the gravita- 
tional constant, M ns the mass of the neutron star and k 
the electron scattering opacity. Then the distance may be 
estimated as: 



d _ f LEdd,non_rel 
.max. 



1/2 



(B.l) 



where £ is the anisotropy factor (for isotropic radiation 
£=1). However, when the photosphere remains rather 
close to the neutron star, gravitational redshift effects be- 
come important (due to relativistic time dilation, see e.g. 
Lewin et al. 1993). The distance derived using Eq. B.l 
is in fact only valid for values of the photosphere radius, 
R»R ns (see below). For R<;R ns Eq. B.l in principle only 
provides an upper limit on d. 

The (relativistic) Eddington luminosity for an observer 
near the neutron star is given by: 



For an observer at Earth, Eq. B.2 becomes, taking into 
gravitational redshift effects: 



L 



Edd. 



L Edd [l-2GM ns /(Rc 2 )] = L Edd [l-R s /R],(B.3) 



where R g is the Schwarzschild radius (R>1.5R 9 ; see Lewin 
et al. 1993). Also the effective black-body temperature 
measured at Earth, T e ff jCO , is affected by the gravitational 
potential: 



T e ff,oc — T c ff [1 — Rg/R] 1 / 2 , 



(B.4) 



where T e ff is the effective black-body temperature near 
the neutron star. Assuming that the source radiates as a 
black-body, we know that: 



L Edd = 47rR 2 <TT e 4 ff . 

For L E dd,oo we also have: 

LEdd.oo = 47T(Z 2 £F E dd,oc, 



(B.5) 



(B.6) 



where F E dd,oo is the Eddington flux received at Earth. 

Eqs. B.3 and B.4 are strictly valid only for non-rotating 
stars, but it is still approximately correct for stars with 
rotation periods of a few milliseconds. F E dd,oo is set equal 
to the observed net-burst fluxes, Fbb,nct, during the radius 
expansion/contraction phase. As noticed in Sect. 5.2, it 
is not easy to infer T ff !00 from our observed values of 
Tbb- We therefore fixed the hardening factor, Tbb/T ff !00 , 
to different values in order to assess its effect. We also 
assume certain values for M ns and £, and are then left 
with five equations (B.2-B.6) and five unknown variables 
(L E dd, L E dd,oc, T c ff, R and d); this enables us to derive a 
distance estimate. First, one can eliminate L E dd and T e ff 
by combining Eqs. B.2, B.4 and B.5, and one gets the 
following equality (for R>1.5R g ): 



cGM n 



R 2 



R 



R — Rc 



(B.7) 



L Edd = (4^cGM ns / K )[l - 2GM ns /(Rc 2 )]- 1 /2. 



(B.2) 



Since this equation can not be solved easily in an analytic 
way, we solved for R numerically. Using R we can then 
determine L E dd,oo (Eq. B.3), which in turn can be used to 
determine the distance, d (Eq. B.6). 

In the low-temperature limit the electron scattering 
opacity, k, is given by k — 0.2(1 + A)cm 2 g _1 , where 
X is the hydrogen fraction (by mass) of the photo- 
spheric matter (see e.g. Lewin et al. 1993). Note that 
for cosmic compositions A=0.73. At very high temper- 
atures (i.e. probably near the peak of bursts), however, 
the scattering electrons become relativistic and k may 
instead be approximated by (in the low-density limit) 
k = k [1 + (kT/39.2 keV) - 86 ]- 1 , where k is same as 
for the low-temperature limit (Paczyhski 1983). It is not 
entirely clear how kT relates to our observed values of 
kTbb, but one may assume kT=kTbb (e.g. van Paradijs et 
al. 1990). We use the high temperature electron scatter- 
ing opacity in our calculations. The estimated distances 
differ by a factor of ~1.05 with respect to the estimates 
when taking the electron scattering opacity in the low- 
temperature limit (see also Table B.2). 
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Table B.l. Distance (d in kpc) and radius at touch down from radius expansion bursts' 







0.25 s spectra 




16 s spectra 


Rtd'' 


burst 


# c 


j d 
, - t non_rcl 




# 


^non_rcl 




(km) 


b4 


84 


11.4±0.3 


10.8±0.3 


4 


11.6±0.3 


10.9±0.3 


18.3 


1)6 


83 


12.1±0.5 


11.4±0.5 


6 


12.2±0.1 


ll.lzbO.l 


12.8 


b7 


94 


12.9±0.6 


12.0±0.5 


2 


12.7±0.2 


11.8±0.1 


16.7 


1)8 


93 


12.3±0.6 


11.5±0.5 


2 


12.1±0.1 


11.2±0.1 


14.4 


1)9 


82 


12.4±0.5 


11.5±0.5 


1 


12.1 


11.1 


14.9 


Averac 




12.2±0.2 


11.4±0.2 




12.1±0.1 


11.2±0.1 


15.4±1.0 



a For the parameters used we refer to the text. 

6 Black-body radius R at touch down for a local observer (see text). 
c Number of spectra in the expansion/contraction phase. 
d Not taking into account gravitational redshift effects. 
e Taking into account gravitational redshift effects. 



The above described procedure for estimating the dis- 
tance was followed using the values for Fbb,net and Tbb for 
each of the spectral fits to the 0.25 s and 16 s spectra in 
the radius expansion/contraction phase. When we assume 
no spectral hardening, i.e. T ff iOO =Tbb, and 'standard' 
parameters (M ns =lAM.Q, X=0.73) no solutions can be 
found for those spectra having kTbb<d-75 keV, i.e. only at 
the largest expansion radii do we have solutions (see also 
Table B.2)^|. The highest observed values of kTbb during 
the radius expansion/contraction phase give the strongest 
constraints on the spectral hardening factor. This is at 
the moment of 'touch-down', i.e. at the vertex of the hor- 
izontal track and diagonal track in the flux versus tem- 
perature diagram. This locus is well defined for our long 
bursts observed in GX17+2. Using the loci of the radius 
expansion bursts b4, b6-b9, we find that the minimum 
value of the (constant) hardening factor is between 1.34 
(burst b4) and 1.52 (burst b6). The distance and radius R 
as inferred at the loci both increase for increasing spectral 
hardening (e.g. for burst b6, the distance increases by a 
factor ~1.1, while R increases by a factor of ~2 when the 
spectral hardening factor increases from 1.6 to 2). 

In Table B.l the results for 'standard' burst parame- 
ters can be seen, i.e. M DS — 1.4 Mq, cosmic composition of 
the photospheric matter, isotropic radiation (£=1), as well 
as taking Tbb/T e ff l0 o=l-7 and using the electron scatter- 
ing opacity appropriate for high temperatures. We show 
the distance as derived by including gravitational redshift 
effects (droi) and not including these effects (d n on_rci, he. 
using Eq. B.l). The errors are set equal to the observed 
variances in the derived distances. 

By relaxing the assumptions of standard burst param- 
eters (i.e. varying M ns , X and £ individually within rea- 
sonable limits, while keeping the others at the 'standard' 
values) one can get an estimate of the systematic uncer- 
tainties involved (see e.g. also Kuulkers & van der Klis 
2000). The results are shown in Table B.2 for burst b6, 

5 Note that the gravitational redshift corrected distance es- 
timates by Smale (1998: CygX-2), Kuulkers & van Klis (2000: 
GX3+1) and Kaptein et al. (2000: 1RXS J171824. 2-402934) 
are derived using T e ff,oo=Tbb and R=Rbb, thus only using 
Eqs. B.2, B.3 and B.6. 



which is the longest burst of our sample. We show the 
results assuming a neutron star mass of M ns = 2 Mq , 
in the range required in models for explaining the kHz 
QPO in neutron star LMXBs (see e.g. Stella et al. 1999; 
Stella 2001; Lamb & Miller 2001), and close to the dy- 
namical neutron mass estimate of CygX-2 using opti- 
cal spectroscopy/photometry, M ns = 1.78±0.23 (Orosz & 
Kuulkers 1999). Note that the distance derived is rather 
high (Table B.2). We also made the assumption that the 
photospheric composition is hydrogen-poor (i.e. X=0), 
as has been argued to be the case during radius expan- 
sion (see Sugimoto et al. 1984). However, the contrac- 
tion phase of long bursts are typical of unstable mixed 
hydrogen/helium burning (see Sect. 5.4), which suggests 
that this is not a reasonable assumption. Also, the per- 
sistent mass accretion rate is inferred to be high, mak- 
ing it plausible that the photosphere will be continu- 
ously supplied with new hydrogen matter. Moreover, the 
distances derived are rather high (see Table B.2). We 
changed the anisotropy factor between the reasonable val- 
ues 0.5 ^ £ ^ 2 (see van Paradijs & Lewin 1987; Lewin 
et al. 1993). Finally, we compared the distance values de- 
rived using values for k valid for low temperatures. The 
results shown in Table B.2 suggests that the systematic 
uncertainties on the derived distance are up to ^30%. 
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Table B.2. Distance (d in kpc) and radius at touch down for burst b6 using different assumptions 





0.25 s spectra 


16 s spectra 


Htd 


parameter 


J c 


d IC \ 


^non_rcl 




(km) 


standard e 


12.1±0.5 


11.4±0.5 


12.2±0.1 


ll.lzbO.l 


12.8 


M ns =2.OM 


14.4±0.6 


13.5±0.6 


14.6±0.1 


12.7±0.2 


13.0 


X=0 


15.9±0.7 


15.3±0.6 


16.1±0.1 


15.1±0.1 


18.9 


low-T K 


11.6±0.5 


11.0±0.5 


11.7±0.1 


10.6±0.1 


11.8 




8.5±0.3 


8.1±0.3 


8.6±0.1 


7.9±0.1 


12.8 


£=0.5 


17.1±0.7 


16.2±0.7 


17.3±0.1 


15.8±0.1 


12.8 


Tbb/T fr,oo=2 


12.1±0.5 


11.7±0.5 


12.2±0.1 


11.6=1=0.1 


19.9 


Tbb/T ff >00 = l^ 


12.1±0.5 (83) 


10.5±0.7 (9) 


12.2±0.1 (6) 







a See text for details. 

6 Black-body radius R at touch down for a local observer (see text). 

c Not taking into account gravitational redshift effects. 

d Taking into account gravitational redshift effects. 

e M ns =1.4M , X=0.73, T bb /T cff , 00 =1.7,high-T k, £=1. 

* Between brackets the number of spectra taken into account are given. 



